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Chapter  1 


Summary 


In  the  years  to  come,  the  size  and  cost  of  cryo-coolers  will  get  smaller 
and  the  demand  for  a  VLSI  high-temperature  superconducting  (HTS)  Josephson 
junction  technology  will  increase.  One  possible  candidate  to  fill  this  need  is  the 
“ion-damage”  Josephson  junction.  These  junctions  are  fabricated  by  using  ion 
bombardment  to  create  localized  narrow  regions  of  defects  in  the  plane  of  a  thin 
film  superconductor.  These  regions  have  a  superconducting  transition  temperature 
lower  than  that  of  the  bulk  film  and  act  as  non-hysteretic  Josephson  junctions.  The 
advantage  of  these  junctions  over  other  technologies  is  that  they  have  no  interfaces 
between  different  materials,  and  can  be  placed  over  10  times  closer  to  each  other 
in  comparison  to  competing  techniques. 

Individual  junctions  and  multi-junction  serial  arrays  were  fabricated  and 
characterized.  Current-voltage  characteristics  were  well  described  by  the  restively 
shunted  junction  model  at  temperatures  close  to  the  critical  temperature  of  the 
weak  link.  At  lower  temperatures  the  junctions  became  strongly  coupled  and 
exhibit  flux  flow  characteristics.  It  is  suggested  that  this  is  due  to  the  barrier 
changing  from  non-superconducting  to  superconducting,  which  allows  Abrikosov 
vortices  to  enter  the  junction.  Junction  arrays  have  flat  giant  Shapiro  steps  at 
N  times  the  voltage  predicted  by  the  ac  Josephson  relation,  where  N  is  the  num¬ 
ber  of  junctions  in  the  array.  Modulation  of  the  critical  current  was  measured  to 
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observe  the  dc  Josephson  effect.  Unlike  other  planar  junction  types  these  junc¬ 
tions  do  not  exhibit  flux-focusing  effects  in  magnetic  field  measurements  and  are 
more  comparable  with  classical  sandwich  type  Josephson  junctions.  Low  temper¬ 
ature  (100°C)  anneals  substantially  increased  the  critical  current  and  operating 
temperature  possibly  from  the  recombination  of  vacancy-interstitial  pairs  formed 
during  implantation.  Junctions  were  also  fabricated  from  the  low  temperature  su¬ 
perconductor  MgB2  for  comparisons  with  the  HTS  junctions  and  to  show  that  this 
technology  is  not  limited  to  HTS. 

Overall  it  has  been  shown  that  this  technique  can  be  used  to  reproduce 
junctions  with  uniform  resistances  and  critical  currents.  With  a  few  innovations 
in  materials  and/or  lithography,  ion  damage  Josephson  junctions  may  become 
common  place  in  the  medical,  communications,  and  defense  industries 
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Chapter  2 

Introduction 

2.1  Historical  Background 

In  1911  Kamerlingh  Onnes  discovered  that  electrical  resistance  of  mer¬ 
cury  wires  abruptly  dropped  to  zero  when  they  were  cooled  below  4.2  K  [1],  He 
named  this  phenomenon  superconductivity  and  in  the  years  to  follow  many  other 
materials  were  found  to  exhibit  this  effect.  Besides  zero  electrical  resistance  Meiss¬ 
ner  found  superconductors  are  perfect  diamagnets  and  expel  magnetic  field  [2]. 
Phenomenological  theories  put  forth  by  the  Londons  [3],  and  Ginsburg  and  Lan¬ 
dau  [4],  described  superconductivity  in  terms  of  electrodynamics  and  thermody¬ 
namics.  This  explained  the  observation  of  the  Meissner  effect  and  revealed  that  the 
transition  from  normal  conductor  to  superconductor  is  a  2nd  order  phase  transition 
which  happens  upon  cooling  below  a  temperature  unique  to  the  specific  material 
called  the  critical  temperature,  Tc.  For  most  conventional  superconductors  the 
mechanism  of  superconductivity  is  well  understood  in  the  context  of  a  microscopic 
theory  by  Bardeen,  Cooper,  and  Schrieffer  (BCS  theory)  [5]. 

2.2  Coherent  Phenomena  in  Superconductors 

The  effects  observed  in  superconductors  are  a  result  of  macroscopic  co¬ 
herence  of  the  current  carriers  in  the  material.  The  motion  of  charge  carriers  in 
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normal  matter  can  be  described  by  the  ordinary  Shrodinger  equation: 

ih ^  =  m  (2.i) 

where  %l>  is  the  wave  function  of  a  particle, 

if)  =  \'>p(r,t)\  exp  [i9(r,t)}  (2.2) 

H  is  a  Hamiltonian,  6  is  th  ephase  of  the  wavefunction,  and  h  is  Planck’s  constant. 
In  the  stationary  state  \ip\  is  constant  because  from  quantum  mechanics  \ip\2  is 
proportional  to  the  density  of  particles.  H  can  be  replaced  by  E  and  equation 
(2.1)  takes  the  form, 

h—  =  -E  (2.3) 

Equation  (2.3)  shows  that  the  quantum  behavior  of  a  carrier  is  reduced  to  a  func¬ 

tion  of  its  wavefunction ’s  phase.  In  normal  conductors,  electrons  obey  Fermi-Dirac 
statistics  and  only  two  carriers  can  have  the  same  energy;  therefore  the  rate  of 
phase  change,  ||  can  only  be  the  same  for  two  carriers.  Macroscopic  observables 
are  determined  by  summing  over  a  large  number  of  particles  and  therefore  phase 
information  from  the  individual  carriers  is  lost. 

In  superconductors,  BCS  theory  shows  that  two  electrons  with  opposite 
momentum  and  spin  can  interact  with  lattice  phonons  and  form  a  bound  pair 
called  a  Cooper  pair.  Cooper  pairs  obey  Bose-Einstein  statistics  because  the  sum 
of  two  electrons’  spin  is  an  integer.  Therefore  many  carriers  may  have  the  same 
energy.  As  a  consequence,  many  of  them  at  the  Fermi  surface  can  condense  into 
a  single  collective  state,  the  superconducting  ground  state.  As  a  result  their  || 
are  identical.  This  results  in  macroscopic  quantities  that  depend  on  phase,  i.e. 
current. 

Only  electrons  close  to  the  Fermi  surface  can  participate  in  pairing.  Elec¬ 
trons  with  higher  or  lower  energies  remain  unpaired  and  are  called  quasi-particles 
to  distinguish  them  from  the  paired  electrons  in  the  superconducting  state.  It  is 
energetically  favorable  for  the  pairs’  wavefunctions,  ip  to  phase  lock  forming  a  su¬ 
perfluid  condensate,  which  can  be  described  by  a  single  wave  function  commonly 
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called  the  order  parameter,  4/.  There  is  phase  coherence  between  the  pairs  over  a 
length  scale  called  the  coherence  length,  £.  In  a  volume  of  £3  there  are  typically 
106  pairs  in  this  coherence  volume  for  conventional  superconductors  (e.g.  Nb,  Pb). 
This  coherence  length  is  one  of  the  two  important  length  scales  in  superconductiv¬ 
ity.  The  phase  coherence  leads  to  two  important  phenomena  in  superconductors, 
magnetic  flux  quantization  and  the  Josephson  effect. 

The  phase  of  the  order  parameter  must  be  unique  and  single- valued,  thus 
around  a  closed  loop  it  must  return  to  itself  (modulo  2n).  This  condition  requires 
any  flux  passing  through  the  loop  to  be  quantized  to  the  values, 

$  =  n$0  (2.4) 

where  n  is  an  integer  and  <h0  is  the  flux  quantum, 

=  2.07  x  10-15  Webers 
2e 

2.3  Superconductors  in  Magnetic  Fields 

Superconductors  repel  magnetic  fields  in  the  superconducting  state  as 
expected  for  a  perfect  conductor;  however  Meissner  found  that  field  was  also  ex¬ 
pelled  from  a  superconductor  while  undergoing  its  phase  transition  from  normal 
to  superconducting.  This  is  contrary  to  what  classical  physics  predicts,  for  a  per¬ 
fect  conductor  should  trap  the  field  inside  [6].  In  the  superconducting  case,  the 
magnetic  fields  generate  superconducting  eddy  currents  that  create  an  exactly  op¬ 
posing  field  that  repels  the  flux  making  it  a  perfect  diamagnet.  Heinz  and  Fritz 
London  showed  that  in  a  bulk  superconductor,  there  is  a  length  scale  for  these 
screening  currents  and  over  that  length  the  magnetic  field  penetrates  the  surface 
and  decays  exponentially  [3], 

B  =  B0  exp 

where  B  is  magnetic  field,  x  is  the  distance  from  the  surface  and  Xl  is  the  London 
penetration  depth,  a  material  specific  constant.  This  screening  length  Xl  is  the 


second  important  length  scale  in  superconductivity. 

In  the  presence  of  a  magnetic  field  the  free  energy  of  the  normal  state 
is  lower  and  the  superconductivity  is  quenched.  This  happens  at  a  field  called 
the  thermodynamic  critical  field,  Hc.  This  places  a  limit  on  the  current  that  may 
be  carried  by  a  superconductor  due  to  the  magnetic  fields  induced  by  the  current 
itself  [6].  This  limit  is  called  the  critical  current,  Ic. 

The  ratio  of  the  two  superconducting  length  scales  is  called  the  Ginsburg- 
Landau  parameter  [4], 


This  value  plays  an  important  role  in  how  superconductors  respond  in  magnetic 
fields.  In  typical  classical  clean  elemental  superconductors  (Type  I),  A/,  ~  50  nm 
and  £  ~  300  nm  resulting  in  k  <  1.  In  this  case  all  magnetic  field  is  expelled 
from  the  interior  for  fields  below  a  critical  field  called  Hc i  and  this  is  the  Meissner 
effect.  Fields  in  excess  of  Hc\  quench  superconductivity  and  return  the  material  to 
its  normal  state.  In  superconductors  where  \L  is  larger  then  £  or  /c  >  (Type  II), 
Abrikosov  showed  that  it  is  thermodynamically  favorable  for  the  superconductor  to 
exist  in  a  mixed  normal  and  superconducting  state  [7].  For  fields  in  excess  of  Hc i, 
the  magnetic  field  penetrates  the  superconductor  in  the  form  of  vortices  where  the 
magnetic  field  is  confined.  Vortices  have  a  normal  core  with  a  radius  of  the  order 
of  £,  where  the  field  lines  penetrate.  The  core  is  surrounded  by  superconducting 
shielding  currents  extending  out  for  a  distance  of  approximately  A l-  Due  to  the 
flux  quantization  in  superconductors  (equation  2.4),  each  vortex  contains  only  one 
flux  quantum,  4>0.  If  field  is  increased  further  beyond  Hc\  more  vortices  form 
to  accommodate  the  additional  flux.  An  additional  increase  to  a  field  called  Hc i 
causes  the  cores  to  overlap  and  the  material  reverts  to  its  normal  state  [6].  Vortices 
in  a  superconducting  current  carrying  wire  or  bridge  (current  density  J)  experience 
a  Lorentz  force  from  the  current, 
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This  force  causes  vortices  to  move  transverse  to  the  current.  The  time  rate  of 
change  of  the  flux  results  in  an  electric  field  parallel  to  the  current, 

E  =  ±f 

c  at 

and  this  acts  like  a  resistive  voltage  [6]. 

It  is  energetically  favorable  for  a  vortex  core  to  exist  in  a  region  in  the 
material  that  is  non-superconducting  such  as  a  defect  or  grain  boundary.  This 
acts  like  a  force  opposing  the  Lorentz  force  from  current  flow  and  is  referred  to  as 
flux  pinning.  For  high  current  applications  in  superconductivity,  defects  or  pinning 
centers  are  sometimes  introduced  into  the  material  to  increase  the  critical  current. 

2.4  Josephson  Junctions 

The  most  significant  phenomenon  related  to  the  phase  coherence  of  the 
charge  carriers  in  superconductors  is  its  ability  to  carry  a  current  at  zero  volt¬ 
age,  a  supercurrent.  A  particularly  interesting  manifestation  of  this  phenomenon 
occurs  when  two  superconductors  are  weakly  connected  and  a  supercurrent  flows 
in  a  similar  fashion.  This  phenomenon  is  called  the  Josephson  effect.  If  two  su¬ 
perconductors  are  coupled  through  a  weak  barrier  such  as  an  insulator  or  normal 
metal  (Figure  2.1)  a  supercurrent  will  flow  in  the  absence  of  a  potential  difference 
(voltage)  if  there  is  a  phase  difference  between  the  order  parameters  of  the  two 
superconductors.  The  difference  between  the  phases  of  the  two  superconductors  is 
defined  as  the  Josephson  phase, 

e2)  .  (2.5) 

The  supercurrent  in  the  simplest  case  is  described  by  the  first  Josephson  equation 

[8], 

Is  —  Ic  sin  <p  .  (2.6) 

Where  the  critical  current,  Ic  is  the  maximum  supercurrent  that  the  junction  will 
sustain  and  is  determined  by  junction  geometry,  temperature,  and  the  materials 
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Insulator  or  Normal  Metal 

Figure  2.1:  Schematic  diagram  of  a  Josephson  junction.  The  order  parameters  from  the  superconducting 
electrodes  decay  into  the  barrier  and  couple  the  wavefunctions  together. 


of  the  weak  link  and  electrodes,  amongst  other  factors. 

Furthermore  if  there  is  a  potential  difference  between  the  two  supercon¬ 
ductors  the  phase  difference  changes  with  time  and  is  described  by  the  second 
Josephson  equation  [8], 


fy  _  =  271 r 

dt  h  $0 


(2.7) 


A  simple  derivation  of  this  was  put  forth  by  Likharev  [9].  Writing  equation  (2.3) 
for  each  of  the  superconductors  and  subtracting  yields, 


dip 

hIt  = 


(2.8) 


The  energy  difference  is  between  the  cooper  pairs  of  charge  2e  in  either  side  and 
can  only  exist  in  the  presence  of  an  electrochemical  potential  difference  V.  Thus 
2eV  can  be  substituted  for  the  energy  difference  and  we  are  left  with  the  second 
Josephson  equation.  The  second  Josephson  equation  is  exact  and  only  depends 
on  fundamental  constants.  If  a  dc  voltage  is  applied  across  a  junction  the  phase 
evolves  with  time  at  a  frequency, 

/  =  jV  .  (2.9) 


This  results  in  an  ac  supercurrent  oscillating  between  the  two  superconductors  at 
frequency  /.  The  proportionality  factor  between  voltage  and  frequency  is  referred 
to  as  the  Josephson  constant, 


Kj  =  ^-  —  483.597879  — 


V 


(2.10) 
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This  forms  the  basis  for  the  definition  of  the  volt  and  the  Josephson  voltage  stan¬ 
dard.  It  is  now  a  common  voltage  standard  and  the  way  the  “volt”  is  main¬ 
tained  [10]. 


2.4.1  The  RCSJ  Model 


An  equivalent  circuit  model  of  a  Josephson  junction  was  put  forth  inde¬ 
pendently  by  Stewart  [11]  and  McCumber  [12].  They  modeled  the  junctions  as 
a  resistor,  capacitor,  and  a  current  source  connected  in  parallel  as  in  Figure  2.2. 
This  is  called  the  Resistive  and  Capacitive  Shunted  Junction  model  (RCSJ).  In 
practice  Josephson  junctions  are  typically  current  biased  and  equating  it  to  the 
currents  through  the  three  parallel  elements  yields, 

V  dV 

I  =  —  +  C— +  Icsinip  (2.11) 


The  first  term  is  the  quasiparticle  current  in  the  finite  voltage  regime.  The  second 
term  reflects  the  capacitance  between  the  electrodes.  The  sinusoidal  term  repre¬ 
sents  the  supercurrent  from  the  first  Josephson  equation  (2.6).  Using  equation 
(2.7)  the  voltage  in  the  resistive  and  capacitive  terms  can  be  rewritten  in  terms  of 
the  Josephson  phase. 


/ 


dV 

2e  dt2  + 


h  dip 

—  +  Ie  sm  <p 


2 eR  dt 


(2.12) 


Dividing  by  Ic,  and  introducing  a  new  time  variable, 


2e 

T  I  yl  l  y  -  t  UJ  J  l\nC 

n 


Figure  2.2:  Equivalent  circuit  for  a  Josephson  junction  used  in  the  RCSJ  model. 
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equation  (2.12)  can  be  simplified  to. 

I  d2p  dp 

Tc=Pcw  +  fr+smv  (2'13) 

where  f3c  is  the  Stewart-McCumber  damping  parameter, 

Pc  =  IcRn  RnC  =  0JjRnC 

n 

and  ujj  is  the  Josephson  frequency.  In  the  case  of  non-hysteretic  junctions,  such 
as  the  ones  studied  in  this  thesis,  the  capacitance  is  small  and  the  capacitive  term 
may  be  omitted.  In  this  limit  equation  (2.13)  becomes  a  first  order  differential 
equation  which  can  be  written  as, 


dp 

dt 


2e 

~K 


IcRn 


sin  p 


(2.14) 


In  this  form  it  is  clear  that  for  I  >  Ic  that  the  instantaneous  voltage  is  positive, 
but  varies  sinusoidaly  with  p.  When  sin^  is  positive  the  supercurrent  is  in  the 
direction  of  the  bias  current  and  the  phase  advances  faster.  When  sin  p  is  negative 
the  supercurrent  is  in  the  opposite  direction  of  the  bias  and  the  phase  advances 
slower.  For  currents  just  above  Ic  the  voltage  across  the  junction  is  a  series  of 
pulses  of  frequency, 


/ 


2e 


Vdc 


(2.15) 


where  V*  is  the  time  averaged  dc  voltage.  These  pulses  have  a  maximum  of  2  IcRn 
and  area  time  integral  exactly  equal  to  the  flux  quantum  such  that, 


/ 


h 

Vdt  =  —  =  $0 

2e 


(2.16) 


The  time  average  dc  voltage  may  be  found  be  integrating  equation  2.14, 


V 


0, 


V  =  ICRn\l[T  )  -1, 


for  I  <  Ic 
for  I  >  Ir 


(2.17) 


This  results  in  a  parabolic  current-voltage  characteristic  at  biases  near  /c,  and 
approaches  Ohms  law  at  higher  biases  (Figure  2.3). 


< 

E 


Figure  2.3:  /  —  V  characteristic  for  a  non-hysteretic  Josephson  junction  with  small  capacitance. 

2.4.2  The  Effect  of  an  AC  Bias  on  Josephson  Junctions 

If  a  Josephson  junction  is  biased  with  both  ac  and  dc  simultaneously, 
steps  occur  in  the  supercurrent  (Figure  2.4)  at  the  precise  values, 

^  (2.18) 


Shapiro  Steps  in  Non-Hysteretic  l-V  Characteristic 


-0.6 


Figure  2.4:  Shapiro  steps  in  the  I  —  V  characteristic  for  a  non-hysteretic  Josephson  junction  simultane¬ 
ously  biased  with  both  ac  and  dc  current 
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a 


where  /  is  the  frequency  of  the  ac  bias  and  n  is  an  integer.  These  were  also 
predicted  by  Josephson  in  his  famous  paper  [8].  They  are  called  Shapiro  steps 
named  after  S.  Shapiro  who  was  the  first  to  observe  them  in  a  junction  [13]. 

Integrating  the  second  Josephson  equation  (2.7)  to  determine  <p  an< 
bining  it  with  the  first  Josephson  equation  (2.6)  results  in 

r  rt 


Is(t)  —  7csin 


J  dt  +  ip o 


(2.19) 


Taking  V(t )  to  be  the  sum  of  a  dc  voltage  and  an  ac  voltage  of  frequency  us> 


2eVa, 


Is(t)  =  Ic  ^  (-1)"  Jn~hfL  s™  ~  w*)*  +  T’o] 


(2.20) 


where  Jn  is  a  Bessel  function  of  the  first  kind.  There  is  only  a  dc  component 
to  the  supercurrent  when  nus  =  loj  or  equivalently  V  =  this  creates  steps 
in  the  I  —  V  characteristic  with  half  width  determined  by,  7cJn^^  (Figure  2.5). 
This  derivation  assumes  that  the  junction  is  voltage  biased  which  is  very  difficult 
to  do  in  practice  because  the  impedance  of  a  junction  is  very  small  compared  to 
typical  sources.  For  the  case  of  a  current  biased  junction  an  analytical  solution  is 
not  possible  and  numerical  calculations  are  necessary  [14].  Numerical  results  for 


Figure  2.5:  Shapiro  step  amplitude  dependence  as  a  function  of  ac  bias. 
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a  current  biased  junction  resemble  the  voltage  bias  case  and  show  a  Bessel-like 
dependence  of  the  step  amplitudes  [15]. 


2.4.3  The  Josephson  Penetration  Depth 


For  the  following  discussions  it  is  essential  to  clarify  the  variables  to 
represent  the  spatial  junction  parameters  to  be  used  for  the  remainder  of  this 
dissertation.  Figure  2.6  is  a  schematic  diagram  of  a  planar  junction  with  current 
flow  along  the  z-axis.  The  dimension  of  the  barrier  in  the  direction  of  the  current 
will  be  referred  to  as  the  junction  length,  l.  The  dimension  of  the  bridge  along  the 
y-axis  will  be  referred  to  as  the  junction  width,  w  and  the  thickness  of  the  film 
in  the  z-direction  will  be  referred  to  as  thickness,  t.  In  the  absence  of  a  magnetic 
field  the  supercurrent  in  a  junction  is  distributed  evenly  across  its  length  provided 
that  its  length  is  smaller  then  a  quantity  called  the  Josephson  penetration  depth 


A j , 


Aj  — 


4>n 


(2.21) 


2ttLqJc' 

where  Jc  and  Lq  is  the  critical  current  and  inductance,  per  unit  length  of  the  junc¬ 
tion  respectively.  The  Josephson  penetration  depth  is  a  screening  length  similar 
to  the  London  penetration  depth  in  bulk  superconductors.  For  classical  sandwich 
type  junctions  it  can  be  shown  [14]  that  the  Josephson  penetration  depth  takes 

the  familiar  form,  _ 

A  J  =  —  (2-22) 


47r2y0Jc(2AL  +  0  ’ 


/ 

I  I  Superconductor 
■  Insulator  or  Normal  Metal 


X 


Figure  2.6:  3-dimensional  representation  of  a  Josephson  junction.  The  current  flow  is  along  the  x-axis, 
and  the  width  of  the  barrier  lies  along  the  y-axis. 
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where  jc  is  the  current  density,  Xi  is  the  London  penetration  depth  of  the  material, 
and  l  is  the  length  of  the  junction  barrier.  This  relation  is  not  applicable  to 
junctions  with  a  planar  geometry  because  the  inductance  is  not  the  same  as  in  the 
classical  sandwich  geometry  [16] .  Tolpygo  and  Gurvitch  [16]  derived  an  expression 
for  planar  junctions  by  making  the  assumption  that  the  inductance  was  dominated 
by  the  kinetic  inductance  of  the  film.  In  the  thin  film  limit  where  the  thickness  of 
the  film  t  is  less  than  the  London  penetration  depth  the  kinetic  inductance  can  be 
taken  as  L  «  Using  this  inductance  results  in  the  expression, 

Wfei  <2’23) 

which  differs  from  the  classical  formula  by  about  a  factor  of  \/w/\l  which  implies 
that  the  Josephson  penetration  depth  in  planar  junctions  is  larger. 

The  Josephson  penetration  depth  plays  an  important  role  in  how  cur¬ 
rent  is  distributed  in  a  junction.  Consider  the  junction  depicted  in  Figure  2.7(a). 
The  current  in  the  electrodes  is  carried  along  the  edges  of  the  bridge  extending  a 


I  I  Superconductor 
■  Insulator  or  Normal  Metal 


Current 


Figure  2.7:  Two  dimensional  cross-section  showing  the  current  density  in  a  long  (a)  and  short  (b) 
Josephson  junction.  Current  in  the  electrodes  is  carried  on  the  edges  and  penetrates  a  distance  of  Al 
into  the  interior.  Inside  of  the  barrier  this  screening  length  is  A j. 
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distance  of  A l  into  its  interior.  Upon  reaching  the  junction’s  barrier  the  current 
penetrates  further  up  to  a  distance  of  A j,  but  is  still  screened  from  the  interior. 
If  the  width  of  the  junction  is  comparable  to  A  j  as  in  Figure  2.7(b),  the  current 
distributes  itself  evenly  across  the  junction’s  length.  Junctions  of  this  type  are 
referred  to  as  “short”  junctions  and  in  the  former  case  of  Figure  2.7(a)  they  are 
called  “long”  junctions.  The  junctions  in  this  thesis  work  are  considered  as  short 
junctions  which  simplifies  the  following  discussion  of  the  effect  of  magnetic  field 
on  the  supercurrent. 


2.4.4  The  Effect  of  Magnetic  Field  on  the  Supercurrent 


The  phase  of  a  “short”  Josephson  junction  is  the  same  along  the  width 
of  its  barrier,  except  for  the  case  when  the  junction  is  in  a  magnetic  field.  Flux 
threading  the  junction  causes  the  phase  to  vary  spatially  which  results  in  a  redis¬ 
tribution  of  the  supercurrent.  Consider  the  junction  depicted  in  Figure  2.6  with  a 
field  Hz  applied  along  the  2-axis.  It  can  be  shown  that  the  phase  varies  along  the 
width  of  the  junction  as, 


dip  —2  el 

i  -  — ■ <2-24> 

/ 

where  l  is  equal  to  the  length  of  the  region  penetrated  by  the  field.  This  consists 
of  the  length  of  the  junction  barrier  and  includes  the  distance  the  field  penetrates 
into  the  electrodes  on  either  side  so  that  l  =  l  +  2Al-  For  constant  field  this  can 
be  integrated  and  substituted  into  the  first  Josephson  equation(2.6)  to  obtain  the 
supercurrent  density  along  the  junction, 


Js  =  Jc  sin 
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(2.25) 


where  ip0  is  an  integration  constant. 

The  critical  current  can  be  found  by  integrating  equation  (2.25)  and 
choosing  ip0  so  that  the  current  is  maximized.  In  the  case  of  a  rectangular  shaped 


junction  with  barrier  width  ukbms  results  in, 
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(2.26) 


This  can  be  written  in  terms 
the  flux  quantum  <J>0  =  f ,  . 


(0 

O 


he  flux  through  the  junction,  $  —  lwfj,oHz  and 


(2.27) 


This  function,  plotted  in  Figure  2.8,  is  the  absolute  value  of  a  sine  function  as  in 
a  single  slit  interference  in  optics,  and  is  a  Fraunhofer  pattern  with  a  periodicity 
of  a  flux  quantum  in  the  junction.  It  is  commonly  used  as  a  test  for  the  Josephson 
effect  and  a  measure  of  the  uniformity  of  the  current  distribution  in  a  junction. 
Figure  2.9  shows  the  distribution  of  current  in  the  junction  for  different  values  of 
flux.  Minima  in  the  sine  function  occur  for  integer  values  of  •Fo,  because  there  is 
an  equal  amount  of  current  in  both  directions.  The  maxima  occur  at  half  integer 
values  and  the  maximum  amplitude  is  reduced  as  the  oscillations  become  smaller. 


Dependence  of  Critical  Current  on  Magnetic  Field 


Figure  2.8:  Critical  current  as  a  function  of  magnetic  flux  for  a  rectangular  junction. 
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Figure  2.9:  Spatial  distribution  of  the  supercurrent  in  a  Josephson  junction  for  different  values  of  magnetic 
flux.  One  can  sec  that  at  the  minima  in  the  sine  function  in  Figure  2.8  (integer  values  of  $o)  that  there 
is  an  equal  amount  of  current  in  both  directions.  For  half  integer  values  of  <f>o  the  net-current  flow  is  in 
one  direction  but  its  magnitude  is  reduced  as  more  oscillations  traverse  the  junction  [20]. 


2.4.5  Proximity  Coupling 

When  a  superconductor  is  in  contact  with  a  normal  metal,  cooper  pairs 
from  the  superconductor  diffuse  into  the  normal  metal  and  quasiparticles  from  the 
normal  metal  diffuse  into  the  superconductor.  They  can  extend  for  hundreds  of 
nanometers  in  either  direction  depending  on  the  coherence  length  in  that  material. 
This  has  the  effect  of  reducing  the  Tc  of  the  superconductor  and  making  the  normal 
metal  superconducting.  This  is  called  the  proximity  effect  and  plays  an  important 
role  in  metallic  barrier  junctions.  These  types  of  junctions  can  be  of  two  types. 
The  first  will  be  referred  to  as  SNS  which  stands  for  superconductor- normal  metal 
superconductor.  For  this  junction  the  barrier  is  a  “true  normal”  metal  which  does 
not  undergo  a  superconducting  transition  at  any  temperature.  SN'S  will  be  used  to 
represent  a  junction  that  has  a  normal  metal  barrier  which  has  a  superconducting 
transition  below  the  temperature  of  operation. 

The  pair  amplitude  for  both  cases,  F(x )  decays  into  the  normal  region 
exponentially  in  the  form  of, 

F(x)  oc  e-fc|x|  (2.28) 

where  k  is  determined  by  the  properties  of  the  material  and  temperature.  For 


“clean”  materials, (by  “clean”  it  is  meant  that  the  electronic  mean  free  path  is 
greater  than  the  coherence  length) 


i-l  _  hyfn 

2kbT 


(2.29) 


where  Vfn  is  the  Fermi  velocity  of  the  normal  metal.  If  the  material  is  in  the  dirty 
limit  where  the  mean  free  path,  le<n  is  much  smaller  than  the  coherence  length 

(le,n  ^  £n)  j 

k-1  =  U  (2.30) 


for  SNS  junctions  and, 


k-1  =  1  + 


(2.31) 


for  SN'S  junctions  where, 

(2-32) 

is  the  normal  metal  coherence  length,  that  can  be  much  longer  than  that  of  the 
superconductor.  Due  to  the  longer  coherence  length  in  most  metals  the  barriers  of 
SNS  and  SN'S  junctions  can  be  much  wider  then  in  the  case  of  an  SIS  tunnel  junc¬ 
tion  where  the  barrier  is  insulating.  The  temperature  dependence  of  the  maximum 
super  current,  close  to  Tc  for  these  junctions  was  calculated  by  de  Gennes  [17]  to 


IC(T)  =  7c0  I  1  -  - 


kl 

sinh  kl 


(2.33) 


where  Ico  is  the  critical  current  at  0  K  and  l  is  the  width  of  the  barrier. 


2. 4. 5.1  Ion  Damage  SN'S  Junctions 

In  high-Tc  materials  the  critical  temperature  can  be  significantly  reduced 
by  disorder  inflicted  from  ion  or  electron  beam  damage.  This  has  been  used  to 
create  SN'S  junctions  within  the  plane  of  a  single  thin  film  [18,19],  by  selectively 
damaging  (reducing  the  Tc )  of  the  region  intended  to  serve  as  a  barrier.  These 
junctions  have  no  interface  between  the  higher  Tc  and  lower  Tc  materials  and  both 
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regions  have  nearly  identical  properties  with  the  exception  of  Tc.  Clarke  proposed 
a  model  where  the  interfaces  can  move  with  temperature  as  a  form  of  Andreev 
scattering,  i.e  the  width  of  the  barrier  gets  smaller  as  temperature  is  reduced  [20]. 
To  account  for  this  movement  “soft”  boundary  conditions  have  been  used  with  the 
classical  de  Gennes  model,  equation  (2.33)  to  explain  experimentally  measured 
critical  current  temperature  dependence  [21].  Katz  and  Dynes  [21]  introduced  a 
temperature  dependant  junction  width  into  the  de  Gennes  model  which  results  in 
a  critical  current  dependence  of  the  form, 

/c(T)a  (l  — I)'  (2-34) 

The  additional  factor  of  T,  from  the  temperature  dependant  width,  fits  well  with 
experimental  data  (Figure  2.10).  Furthermore  they  also  introduced  a  current  de¬ 
pendant  resistance  because  at  higher  biases  the  interfaces  move  away  from  the 
center  of  the  weak  link  as  a  local  Ic  is  exceeded  for  more  of  the  material.  They 
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Figure  2.10:  Critical  current  dependencies  for  3  ion-damage  Josephson  junctions  of  different  lengths. 
The  temperature  dependant  barrier  length  provides  a  better  fit  then  the  classical  de  Gennes  model  with 
and  without  “soft”  boundary  conditions  [27]. 
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Voltage  (mV) 

Figure  2.11:  I  —  V  characteristic  for  an  ion-damaged  Josephson  junction  fit  with  the  RSJ  model  modified 
with  a  current  dependant  resistance  [27], 

suggested  a  resistance  of  the  form, 

R(I)  =  R0  +  aI  (2.35) 

where  a  is  a  fitting  parameter.  Substituting  this  into  the  RSJ  model  fit  their  data 
well  (Figure  2.11)  and  explained  a  downward  curvature  observed  in  the  I—V  curve 
for  junctions  of  this  type  [21]. 

2. 4. 5. 2  SS'S  Junctions 

Another  type  of  junction  worthy  of  mention  is  referred  to  as  an  SS'S 
junction,  where  the  barrier  material  like  in  the  SN'S  case  is  also  a  reduced  Tc 
superconductor.  The  difference  between  the  two  is  that  the  SS'S  junction  operates 
below  the  Tc  of  the  barrier  whereas  the  SN'S  junction  operates  above  Tc.  The 
implications  of  having  a  superconducting  barrier  will  be  discussed  in  more  detail 
in  chapter  5. 

2.4.6  Series  Arrays  of  Josephson  Junctions 

Series  or  parallel  arrays  of  junctions  are  desirable  for  applications  where 
junctions  are  used  as  oscillators  [22],  voltage  standards  [23]  or  sensitive  detectors.  A 


21 


Josephson  oscillator  involves  applying  a  dc  voltage  across  the  junction  and  coupling 
the  ac  supercurrent  generated  into  an  antenna.  It  is  difficult  to  do  this  in  practice 
with  single  junctions  because  the  single  junction  impedances  are  typically  very 
small.  Utilizing  series  arrays  of  junctions  increases  the  impedance  and  circumvents 
this  difficulty.  Furthermore  the  power  output  from  an  array  increases  as  the  square 
of  the  number  of  junctions  in  the  array.  In  single  junction  oscillators  the  radiation 
line  width  is  proportional  to  the  normal  state  resistance,  A  to  oc  Rn.  In  the  case  of 
an  impedance  matched  array  it  decreases  by  the  square  of  the  number  of  junctions 
in  the  array,  i.e  Aw  a  [14], 

For  voltage  standards,  arrays  are  used  to  increase  the  voltage  available 
[14].  For  a  single  junction  driven  at  20  GHz  the  voltage  of  the  first  Shapiro  step 
occurs  at  40  /iV  which  is  too  small  for  practical  applications.  However  in  an  array 
the  voltage  from  equation  (2.6)  increases  as  the  number  of  junctions  in  the  array. 


nNhf 

2e 


(2.36) 


Where  n  is  the  step  number,  N  is  the  number  of  junctions  in  the  array  and  /  is 
the  drive  frequency,  so  a  1000  junction  array  at  20  Ghz  yields  40  mV,  a  usable 
value. 


2.4.6. 1  The  Josephson  Arbitrary  Waveform  Generator 

Most  Josephson  voltage  standards  are  dc  standards.  A  junction  array  is 
biased  with  a  high  frequency  ac  current  to  create  constant  voltage  steps.  The  num¬ 
ber  of  junctions  needed  to  produce  the  desired  voltage  are  dc  biased  simultaneously 
with  the  ac  drive  signal.  The  dc  current  is  selected  to  occur  at  the  constant  voltage 
step.  The  frequency  of  the  ac  current  and  the  number  of  junctions  biased  deter¬ 
mine  the  voltage  of  the  Shapiro  step.  This  type  of  standard  may  not  be  used  for 
ac  voltages,  especially  for  high  frequency  (>lMHz)  by  simply  switching  between 
the  desired  voltages  to  create  a  waveform  because  the  voltage  across  the  array  is 
undefined  during  the  switching  process  [24].  Northrop  Grumman  and  NIST  co- 
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Figure  2.12:  Voltage-current  characteristic  for  a  junction  array  simultaneously  ac  and  dc  biased  depicting 
optimum  dc  bias  currents,  ±/o  for  selection  of  the  output  voltage  polarity  [30]. 


developed  a  new  type  of  ac  standard  that  involves  closely  spaced  junctions  in  an 
array  as  a  lumped  circuit  element  [25].  This  circuit  is  called  a  Josephson  arbitrary 
waveform  generator.  It  is  operated  by  biasing  the  junction  with  a  sinusoidal  drive 
signal  and  a  2-level  binary  bit  stream.  Waveforms  may  be  stored  in  the  form  of 
binary  code  for  the  playback  of  waveforms  with  quantum  accuracy. 

Figure  2.12  shows  the  V  —  I  characteristic  of  a  junction  biased  simulta¬ 
neously  with  ac  and  dc  current.  The  voltage  (V  =  n^-)  of  the  n  =  ±1  Shapiro 
steps  are  shown  and  are  the  time  averaged  voltage  of  a  set  of  pulses  whose  time 
integral  is  precisely  equal  to  the  flux  quantum  (as  explained  in  section  2.41). 
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Figure  2.13:  Example  of  dc  voltage  synthesis  using  a  pulse  driven  waveform  generator.  The  top  curve 
is  the  ac  bias  drive  signal,  the  middle  is  the  binary  bit  stream  and  the  bottom  shows  the  instantaneous 
voltage  across  the  junction  [30]. 
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Figure  2.14:  2-level  binary  code  used  to  generate  the  sine  wave  S' ,  the  pulses  in  So  represent  the  density 
of  l’s  and  0’s  in  the  bit  stream  [30]. 

The  bias  currents  +/0  and  —I0  are  used  to  select  the  polarity  of  the  voltage  out¬ 
put.  Time  averaged  voltages  between  these  two  values  can  be  generated  by  using 
periodic  sequences  of  bias  pulses.  Figure  2.13  shows  an  example  of  a  6  bit  code. 
The  top  curve  is  the  drive  signal  used  to  ac  bias  the  array.  The  middle  shows  the 
binary  pattern  used  to  select  which  voltage  step  is  output,  and  the  third  trace  is 
the  time  dependant  output  of  the  junction  array.  If  this  6  bit  pattern  is  continu¬ 
ously  repeated  it  results  in  2  extra  positive  polarity  pulses  for  each  cycle  resulting 
in  a  voltage  equal  to  V  —  Similarly  ac  voltages  can  be  generated  using 

complex  bit  patterns.  Figure  2.14  shows  a  two  level  binary  code  used  to  create  a 
sine  wave.  This  code  represents  the  density  of  the  l's  and  0's  in  the  bit  stream.  In 
the  peak  there  is  a  larger  density  of  l's  and  in  the  node  there  are  an  equal  number 
of  l's  and  0's  (i.e.  the  number  of  pulses  of  positive  polarity  equals  that  of  negative 
polarity  resulting  in  a  zero  voltage). 

For  this  device  it  is  essential  that  the  Josephson  array  appears  as  a  lumped 
element  to  the  drive  signal  and  the  bit  stream  so  that  each  junction  in  the  array 
sees  the  same  power.  By  “lumped  element”  it  is  meant  that  the  device  is  smaller 
then  \  wavelength  of  the  microwave.  For  a  1  volt  standard  driven  at  18  GHz  Benz 
has  determined  the  optimum  junction  parameters  for  this  application  [23].  It  would 
require  13,500  junctions  with  37/zF  IcRn  products  and  an  inter-junction  spacing 
of  120  nm.  The  junctions  to  be  described  in  this  thesis  work  were  motivated  by 
this  application. 
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Figure  2.15:  Circuit  diagram  for  a  dc  SQUID. 


2.4.7  The  DC  SQUID 

The  dc  SQUID  [26]  (Figure  2.15)  consists  of  two  Josephson  junctions 
connected  in  parallel  in  a  superconducting  loop  with  inductance  L.  Biasing  the 
SQUID  with  a  constant  current  IB,  causes  the  voltage  V  across  the  SQUID  to 
oscillate  with  a  period  /0,  as  external  magnetic  flux  $  is  changed  (figure  2.16). 
Small  changes  in  magnetic  flux,  ($  <  <F0)  can  be  measured  by  selecting  a  bias 
current  that  maximizes  the  amplitude  of  the  voltage  modulation  and  sets  the 
external  flux  at, 


i 


Figure  2.16:  Modulation  of  the  current-voltage  characteristic  for  a  dc  SQUID  the  trace  on  the  left  shows 
how  Ic  changes  with  applied  flux  and  the  trace  underneath  shows  the  voltage  modulation  or  “triangle 
wave”  resulting  [20]. 
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(2n  +  l)^,  (71  =  0,1,2,...) 


so  that  the  flux-to-voltage  transfer  coefficient, 

\dV 


U*  = 


9$ 


is  maximized.  Thus  the  SQUID  produces  a  maximum  output  voltage  signal, 


3V  =  149$ 


in  response  to  a  small  flux  signal  9$. 

The  sensitivity  of  the  SQUID  to  detect  small  magnetic  fields  is  ultimately 
limited  by  intrinsic  white  noise  associated  with  the  resistances  of  the  junctions. 
Sensitivity  can  be  improved  by  connecting  many  SQUID  loops  together  in  a  se¬ 
ries  array  [14].  The  voltage  signal  increases  by  a  factor  of  N  while  the  noise  only 
increases  by,  \/N  where  N  is  the  number  of  SQUID  loops  in  the  array.  Nonuni¬ 
form  junction  parameters  and/or  loop  size  variations  degrade  the  performance  of 
SQUID  arrays.  That  has  made  them  difficult  to  realize  in  high-Tc  superconducting 
materials. 


2.5  Superconducting  Materials 

The  first  materials  discovered  to  exhibit  superconductivity  were  elemental 
superconductors  such  as  Hg,  Pb,  A1  and  Nb.  Clean  elemental  superconductors 
are  type  I  superconductors  with  the  exception  of  Nb.  They  have  relatively  long 
coherence  lengths  and  low  critical  temperatures  with  the  highest  being  Nb  at  9.25 
K.  Niobium  is  commonly  used  for  low-temperature  Josephson  junctions  for  SQUIDs 
and  digital  Josephson  applications.  Historically,  soon  to  follow  were  numerous 
binary  compounds  that  exhibited  superconductivity.  Some  of  these  being  type  II 
making  high  field  magnets  possible.  Some  of  the  more  notable  ones  are,  NbN, 
Nb3Sn,  and  Nb3Ge.  The  latter  two  belong  to  a  group  of  compounds  called  the 
A15  superconductors  named  after  their  crystal  structure.  This  pushed  the  Tc  up 
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to  23.2  I<  (Nb3Ge)  a  record  which  held  for  binary  compound  superconductors  until 
the  discovery  of  MgB2  in  2001  with  a  Tc  of  39  K.  All  of  these  superconductors  are 
described  well  by  BCS  theory. 

In  1986  Bednorz  and  Muller  discovered  superconductivity  in  a  layered 
ceramic  material  I^-zBazCuC^  [27].  Its  Tc  of  30  K  was  much  higher  then  any 
other  superconductors  at  the  time  and  this  discovery  sparked  a  large  research  effort 
in  the  new  materials.  In  1987  a  major  breakthrough  came  with  the  discovery 
of  YBa2Cu307_i  [28].  Its  Tc  of  92  K  was  the  first  superconductor  to  have  Tc 
above  the  boiling  point  of  liquid  nitrogen  77  K.  This  was  a  major  breakthrough  as 
this  potentially  simplified  cooling  requirements  and  reduced  cost  for  a  variety  of 
applications. 

2.5.1  YBa2Cu307_s 

YBa2Cu307_,5  (YBCO)  belongs  to  a  family  of  oxygen  deficient  pervoskite 
layered  materials,  R-Ba2Cu307_<5,  where  R  is  one  of  the  Rare  earth  elements,  Y, 
La,  Pr,  Nd,  Sm,  Eu,  Gd,  Tb,  Dy,  Ho,  Er,  Tm,  Yb,  and  Lu.  YBa2Cu307_i  is  the 
most  studied  and  widely  used  of  these  compounds  and  the  primary  material  used 
in  this  thesis.  The  nature  of  superconductivity  in  these  materials  is  not  explained 
by  BCS  and  has  eluded  physicists  for  almost  2  decades,  however  empirical  obser¬ 
vations  have  provided  some  clues.  Figure  2.17  shows  a  single  unit  cell  of  YBCO. 
It  may  exist  in  two  different  structures  depending  on  the  reduced  Oxygen  content 
(5).  With  <5=1  the  0(1)  sites  in  the  Cu(l)  layer  are  vacant  and  the  structure  is 
tetragonal.  In  this  form  YBCO  is  insulating  and  does  not  superconduct.  Increas¬ 
ing  the  oxygen  content  adds  atoms  at  the  0(1)  sites.  When  5  <  0.65  oxygen  atoms 
at  the  0(1)  sites  form  Cu-0  chains  along  the  6-axis.  This  elongates  the  6-axis  and 
changes  the  structure  to  orthorhombic,  with  o=3.82  A,  6=3.89  A,  and  c=11.68 
A.  Optimum  superconducting  properties  occur  when  <5  ~  0.07  and  all  of  the  0(1) 
sites  are  occupied  with  the  exception  of  a  few  vacancies.  Substitution  experiments 
of  atoms  at  the  Cu  and  rare  earth  sites  provide  evidence  that  conduction  occurs  in 
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Figure  2.17:  Orthorhombic  unit  cell  for  YBa2Cu307-j,  showing  the  positions  of  the  different  atoms. 
The  cell  is  usually  described  as  three  oxygen  deficient  perovskite  cells  stacked  on  top  of  each  other. 


the  Cu(2)0  planes  and  the  Cu(l)0(l)  chains  act  like  a  charge  reservoir,  which  pro¬ 
vides  carriers  to  the  CuO  planes.  With  conductivity  confined  to  a- b  planes  a  large 
anisotropy  in  transport  properties  results.  Along  the  c-axis,  normal  conductivity 
is  10  times  smaller  then  in  the  a—b  plane.  For  other  cuprates  in  the  same  general 
class,  the  anisotropy  is  even  greater  and  inter  plane  transport  is  highly  restricted. 
Furthermore  the  superconducting  length  scales  show  similar  anisotropy,  the  pene¬ 
tration  depth  ( Xab  «  150  nm,  Ac  «  800  nm)  and  the  coherence  length,  (£ab  «  2  nm, 
£c  «  0.4  nm).  Although  the  coherence  length  in  the  a  —  b  plane  is  5  times  bigger 
than  that  along  the  c-axis  it  is  quite  small  compared  to  classic  superconductors 
such  as  niobium,  (f  ss  40  nm).  This  means  that  the  superconducting  state  is  more 
susceptible  to  local  disruptions  from  interfaces  or  defects  of  the  order  of  a  single 
unit  cell.  This  complicates  device  fabrication.  Another  difficulty  in  working  with 
this  material  is  that  it  readily  forms  twin  boundaries  and  is  easily  degraded  by 
water. 

2.5.2  MgB2 

The  structure  of  MgB2  has  a  hexagonal  unit  cell,  with  a  =  b  =  3.086  A 
and  c  =  3.524  A.  HCP  layers  of  Magnesium  alternate  with  graphitic  Boron  layers 
(Figure  5.2). 
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In  its  normal  state  it  is  metallic  and  has  a  room  temperature  resistivity  of 
about  5  pflcm.  The  Tc  of  bulk  samples  is  39  K  [29]  and  has  been  reported  as  high 
as  41  K  in  thin  films  [30].  The  penetration  depth  is  around  100  nm  [30]  and  like 
YBa2Cu307_^  the  coherence  length  is  anisotropic.  For  thin  films  the  coherence 
lengths  were  measured  to  be,  £a{,  «  5  nm  and  £c  «  3  nm  [30].  The  higher  Tc 
makes  MgB2  a  promising  candidate  to  replace  niobium  in  commercial  junction 
technologies.  A  recent  survey  of  235  cryocoolers  [31]  reveals  many  potential  benefits 
from  MgB2  devices.  For  a  300  mW  cooler,  the  input  power  to  the  cooler  can  be  up 
to  100  times  less  for  an  MgB2  device  operating  at  37  K  compared  to  Niobium  at 
4.2  K.  Furthermore,  the  cooler  mass  decreases  by  10  times  and  the  cost  is  halved. 


Figure  2.18:  Hexagonal  structure  of  MgB2. 

2.5.3  Ion  Damage  in  Superconductors 

The  superconducting  state  in  some  materials  is  sensitive  to  disorder  which 
can  completely  destroy  superconductivity.  Many  experiments  in  the  A15  com¬ 
pounds  [32]  and  high-Tc  cuprates  [33]  involving  disorder  induced  by  ion  bombard¬ 
ment  have  provided  insight  into  the  mechanisms  involved.  If  the  superconducting 
pair  wave  function  is  written  as, 

9  =  5e'v 

where  S  is  the  pair  amplitude  and  <p  is  the  phase,  it  is  evident  that  there  are 
two  ways  to  destroy  superconductivity;  either  by  reducing  the  pair  amplitude  or 
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destroying  the  phase  coherence.  In  A15  compounds  it  is  thought  that  ion  damage 
creates  defects  that  alter  the  pair  density  whereas  in  cuprates  the  defects  affect  both 
the  phase  coherence  and  the  pair  amplitude.  In  the  A15’s  ion  damage  decreases 
the  critical  temperature  until  it  saturates  when  the  damage  level  decreases  the 
mean  free  path  to  inter-atomic  distances  [33].  Electron  defect  scattering  lowers 
the  density  of  states  and  reduces  the  pair  amplitude  which  results  in  the  lowering 
of  Tc.  In  the  cuprates  there  is  strong  evidence  that  the  mechanism  of  destruction 
of  the  superconducting  state  is  mainly  due  to  defects  in  the  oxygen  ordering  which 
create  non-superconducting  grains  that  disturb  the  phase  coherence  and  the  pair 
density  [33,34], 

2. 5. 3.1  Ion  Damage  in  YBaaCusO?-^ 

Many  studies  of  ion  damage  in  YBa2Cu307_5  have  been  carried  out  over 
the  past  20  years  [33,34].  It  has  been  shown  that  as  dose  increases  Tc  is  reduced 
and  resistivity  is  increased.  For  high  doses,  films  become  highly  disordered  and 
insulating  and  superconductivity  is  destroyed.  Figure  2.19  shows  an  example  of 
the  metal  insulator  transition  of  YBa2Cu307_<5  irradiated  with  1  MeV  Ne+.  In  this 
thesis  non-hysteretic  junctions  are  created  by  the  selective  masking  of  a  bridge  of 
the  superconductor  and  ion  damaging  the  material  with  a  dose  that  lowers  the 
film’s  Tc  but  leaves  it  metallic. 
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Figure  2.19:  Resistivity  as  a  function  of  temperature  for  a  single  YBa2Cu307-i  film  after  bombardment 
with  1  MeV  Ne+  ions  at  fluences  of  (0,  0.1,  2.5,  4.0,  10.0,  15.0,  20.0,  and  22.0)xl013  ions/cm2  [40]. 


2. 5. 3. 2  Ion  Damage  in  MgB2 

A  recent  study  of  ion  damage  in  MgB2  reveals  that  the  mechanism  of 
disorder  is  similar  to  that  of  A15  materials  [35].  The  pair  amplitude  is  decreased 
by  inter-grain  scattering.  For  high  doses  Tc  is  reduced  to  about  7  K  and  there 
is  no  insulator  transition  as  in  the  case  for  YBa2Cu307_,5.  Figure  2.20  shows  the 
results  of  a  sample  bombarded  with  2  MeV  alpha  particles.  Josephson  junctions 
using  this  material  will  also  be  presented  in  this  thesis. 


Figure  2.20:  Resistivity  as  a  function  of  temperature  for  a  MgB2  film  bombarded  with  various  fluences 
of  2  MeV  alpha  particles  [41]. 


Chapter  3 


Methods,  Assumptions, 
Procedures 

3.1  HTS  Thin  Films 

The  majority  of  the  devices  studied  were  fabricated  from  YBCO  thin 
films.  Initially  films  used  were  grown  in  the  lab  at  UCSD  using  pulsed  laser  ab¬ 
lation.  However  since  this  technique  is  limited  to  producing  small  quantities  of 
samples  per  growth-run,  the  source  material  was  changed  to  commercial  films  de¬ 
posited  on  2”  wafers  by  thermal  co-evaporation.  The  commercial  films  proved 
more  economical  and  more  importantly  allowed  the  fabrication  of  many  devices 
from  more  uniform  material  for  comparison  of  other  process  variables. 

For  this  work,  the  majority  of  the  substrates  used  were  CeC>2  buffered  r- 
plane  sapphire.  This  choice  was  made  to  reduce  costs  and  to  ensure  future  devices 
could  be  easily  incorporated  into  microstrip  transmission  lines  for  RF  devices. 
Sapphire  is  a  hexagonal  structure  with,  a  —  4.759  A,  and  c  =  12.97  A.  When 
cut  along  the  (1102)  plane  (r-plane)  it  has  a  pseudo- rectangular  surface  structure. 
Ce02  is  a  cubic  structure  that  grows  on  the  r-plane  of  sapphire  with  its  principle 
axes  along  the  [1011]  and  [1210]  directions.  The  lattice  parameter  of  Ce02  is  5.41  A. 
This  is  a  close  match  to  the  a  —  b  plane  diagonal  of  YE^CuaCb^  and  encourages 
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growth  of  c-axis  (001)  YBCO  with  the  [110]  direction  along  the  [100]  direction  of 
Ce02. 

The  electrical  transport  properties  of  these  materials  were  determined  by 
the  photolithographic  patterning  of  a  4-point  bridge  and  measuring  the  resistance 
as  a  function  of  temperature.  This  will  be  explained  in  detail  in  later  sections. 
From  this  measurement  resistivity  can  be  calculated  using  the  bridge  geometry. 
The  y-intercept  from  a  linear  fit  of  the  data  above  the  transition  is  how  the  residual 
resistivity  was  determined,  (p0).  This  provides  information  about  the  quality  of 
the  normal  conductivity.  The  width  of  the  superconducting  transition  provides 
a  measure  of  the  distribution  of  the  critical  temperatures  throughout  the  bridge. 
A  wide  transition  can  be  attributed  to  inhomogeneous  material  with  regions  of 
different  Tc.  Furthermore,  the  temperature  dependent  resistivity  is  different  for 
the  Cu-0  chains  along  the  6-axis  ( pt )  then  the  Cu-0  planes  along  the  a-axis  (pa). 
Measurements  in  de-twinned  single  crystals  [36]  reveal  that  po  has  a  T2  dependence 
while  pa  is  linear  with  T.  Furthermore  pa  is  of  the  order  of  300  pfl-cm  at  room 
temperature  with  poa  near  0,  while  pb  is  3  times  less  and  p0&  can  be  as  large  as  100 
pfi-cm.  Most  films  contain  grains  of  both  a  and  b  axis  orientation  with  respect  to 
the  current  flow,  thus  they  have  resistivity  characteristics  that  are  a  mix  between 
these  two  extremes. 

Figure  3.1  shows  a  typical  result  for  a  Theva  YBa2Cu307_,5  film.  The 
data  show  a  mix  between  a  and  b  axis  transport.  The  small  width  of  the  super¬ 
conducting  transition  infers  that  the  material  is  homogeneous  and  doesn’t  contain 
mixed  phases.  Figure  3.2  shows  the  results  from  a  NdBCO  film.  The  non-linear 
behavior  in  the  normal  region  and  high  resistivity  indicate  poor  quality  material 
and  the  large  transition  width  of  5  K  indicates  that  there  are  mixed  phases  with 
variation  of  critical  temperature.  To  rule  out  processing  damage,  measurements 
were  repeated  using  a  simple  4-point  VanderPaaws  resistance  measurement  with¬ 
out  removing  the  in-situ  contact.  The  Tc  and  transition  width  was  identical  to  the 
patterned  bridge. 
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Theva  1 0Onm  YBC0/AI203 


Temperature  (K) 


Figure  3.1:  Resistance  vs.  temperature  measurement  of  a  Theva  co-evaporated  100  nm  thick  YBCO  thin 
film  on  Ce02  buffered  r-planc  sapphire  substrate.  The  Tc  of  this  film  is  86  K  and  the  residual  resistivity 
is  17  fiClcm  The  inset  shows  the  transition  width  to  be  about  0.5  K. 

3.2  Mask  Design 

The  superconducting  films  were  patterned  with  photolithography  to  cre¬ 
ate  contacts,  electrodes,  and  alignment  marks  for  electron  beam  lithography  (EBL). 


Theva  lOOnm  NCIBCO/AI2O3 


Temperature  (K) 


Figure  3.2:  Resistance  vs.  temperature  measurement  of  a  Theva  co-evaporated  100  nm  thick  NdBCO 
thin  film  on  Ce02  buffered  r-plane  sapphire  substrate.  The  Tc  of  this  film  is  71  K  and  the  residual 
resistivity  is  approximately  300  ^flcm  The  inset  shows  the  transition  width  to  be  about  5  K. 
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Figure  3.3  shows  one  of  the  masks  designed  for  the  superconducting  layer.  The 
light  lines  (green)  are  alignment  marks  for  electron  beam  lithography  and  the  dark 
lines  (blue)  are  part  of  the  circuit.  This  particular  pattern  has  16  individual  4 
fim  wide  bridges  for  4-point  measurements.  Two  contact  pads  are  common  for  all 
bridges,  V~  and  I~.  Each  bridge  also  has  2  unique  additional  terminals  referred 
to  as  V+  and  I+ .  One  feature  unique  to  this  mask  is  that  all  of  the  devices  for  the 
chip  are  contained  in  a  90  x  90  gm2  area  in  the  center  of  the  chip.  This  is  the  field 
of  view  that  our  EBL  system  can  write  without  having  to  be  refocused  or  have  the 


substrate  stepped. 


Figure  3.3:  Sixteen  bridge  photolithographic  graph  mask  dark  iines(blue)  represent  the  circuit  areas  and 
light  lines(grecn)  are  EBL  alignment  marks.  16  bridges  are  located  within  a  90  pm  x  90  pm  central 
region  of  the  mask. 
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3.3  Pattern  Photolithography 

Wafers  were  diced  by  Theva  into  5x5  mm2  samples.  They  were  spin- 
coated  with  a  1.8  pm  layer  of  Shipley  1818  positive  photoresist  spun  at  4000  rpm 
for  30  s.  Samples  were  then  soft-baked  on  a  hotplate  at  90°C  for  2  minutes  to 
remove  the  solvent  from  the  resist.  After  baking,  samples  were  first  exposed  for 
1  min  at  12.5  mW/cm2  in  a  JBA  contact  mask  aligner  using  a  solid  3.6  mm  x 
3.6  mm  square  that  covered  the  resist  in  the  center  but  left  the  edges  exposed. 
Samples  were  developed  in  Shipley  MF-319  until  the  resist  on  the  edges  was  gone. 
This  step  was  essential  because  the  resist  on  the  edges  of  the  samples  was  much 
thicker  then  in  the  center,  which  prevented  good  contact  at  the  essential  circuit 
areas  in  the  center  of  the  sample.  The  pattern  was  then  exposed  using  a  chrome- 
on-quartz  mask  with  a  0.25  pm  tolerance  for  18  s  at  12  mW/cm2.  Exposed  samples 
were  developed  for  approximately  18  s  in  Shipley  MF-319  developer.  Samples  were 
agitated  during  development  to  ensure  that  resist  coming  off  the  sample  did  not 
remain  near  the  pattern  which  can  prevent  development  of  some  areas.  Samples 
were  rinsed  in  deionized  water  to  stop  the  development  and  then  blown  dried  with 
compressed  nitrogen.  The  Patterns  were  examined  with  an  optical  microscope  to 
ensure  there  were  no  defects  or  resist  remaining  to  be  developed.  After  inspection 
samples  were  placed  on  a  hot  plate  and  hard-baked  at  120°C  for  5  minutes  to  cross 
link  the  polymers  in  the  resist  and  make  it  more  resistant  to  ion-milling. 

3.4  Argon  Ion  Milling 

Argon  ion  mill  dry  etching  has  become  the  standard  way  to  pattern  HTS 
samples.  It  involves  subjecting  the  sample  to  an  energetic  beam  of  ions  that 
physically  sputters  the  atoms  off  of  the  material  being  etched.  A  discharge  from 
a  tungsten  filament  ionizes  atoms  in  the  surrounding  gases  and  creates  a  plasma. 
The  ions  in  the  plasma  are  accelerated  and  focused  using  two  grids  placed  above 
the  discharge.  A  potential  difference  is  applied  between  the  grids  which  cause  the 
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ions  to  accelerate  from  the  lower  to  the  upper  grid.  Upon  reaching  the  upper  grid 
the  ions  pass  through  the  holes  in  it  and  exit  the  source. 

A  custom  built  cryogenic  ion  milling  system  consisting  of  a  diffusion 
pumped  bell  jar,  a  3.5  cm  Kaufmann  ion  source,  Commonwealth  Scientific  ID 
2500  ion  drive,  and  a  liquid  nitrogen  cooled  sample  holder  was  used.  Samples  were 
then  attached  to  a  copper  plate  with  Dow-Corning  vacuum  grease  and  loaded  into 
the  bell  jar.  The  chamber  and  gas  lines  were  pumped  to  2  x  10~6  Torr,  and  then 
liquid  nitrogen  from  a  pressurized  Dewar  was  forced  through  a  coil  of  copper  tub¬ 
ing  attached  to  the  sample  holder.  After  the  sample  holder  was  cooled  to  77  K, 
argon  gas  was  introduced  into  the  chamber  through  a  variable  leak  valve  and  the 
pressure  was  adjusted  to  be  2.5  x  10-4  Torr.  The  ion  mill  power  supply  was  then 
turned  on  and  ran  in  auto  mode  with  the  acceleration  voltage  set  to  400  V  with 
a  beam  current  of  5  mA.  This  etched  YE^CusCU-i  at  75  A/min  and  gold  at  400 
A/min.  The  cryogenic  cooling  helps  protect  the  superconductor  from  heating  that 
would  result  in  oxygen  loss  and  it  also  creates  steeper  edges.  After  ion  milling,  the 
resist  was  removed  using  acetone  and  ultrasound,  followed  by  a  rinse  in  isopropyl 
alcohol  and  then  blown-dry  with  compressed  nitrogen 

3.5  Lift-off  and  Gold  Removal 

Another  photolithographic  step  followed  ion  milling  to  create  metallic 
contacts.  Two  different  processes  were  followed  depending  on  whether  the  sample 
already  had  a  contact  layer  deposited  in-situ  or  whether  it  was  a  bare  film  requiring 
metallization.  Samples  requiring  contacts  were  coated  with  resist  as  described  in 
section  3.4.  A  mask  of  the  contact  pads  was  aligned  to  the  sample  and  exposed 
for  22  s  at  12  mW/cm2.  After  exposure  they  were  then  soaked  in  toluene  or 
chlorobenzene  for  60  s.  This  slows  down  the  development  of  the  resist  on  the  surface 
and  creates  an  overhang  feature  (Figure  3.4).  This  feature  prevents  contacts  from 
pulling  off  during  resist  lift-off,  by  ensuring  the  metal  is  discontinuous  at  the  edges. 
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Figure  3.4:  (a)Ge  film  evaporated  on  top  of  a  lift-off  over  hang  created  using  a  chlorobenzene  soak.  The 
film  is  discontinuous  at  the  resist  edge  preventing  it  from  being  pulled  off  when  the  resist  is  removed, 
(b)  Cross-section  showing  the  extent  of  the  overhang. 

After  development  the  sample  was  ashed  in  oxygen  plasma  in  a  Technics  500  II 
asher  at  200  W  for  30  s  at  a  pressure  of  200  mTorr.  This  descums  the  surface 
and  improves  contact  adhesion.  After  ashing,  a  layer  of  gold  was  deposited  via  dc 
sputtering  in  a  Denton  Discovery  18  sputtering  system.  Then  the  underlying  resist 
was  removed  using  acetone  and  ultrasound. 

For  samples  already  having  a  contact  layer  deposited  in-situ,  the  metal 
over  the  center  of  the  sample  had  to  be  removed.  This  was  done  by  exposing  and 
developing  a  500  x  100  fx m2  area  over  the  16  bridge  region  (Figure  3.5).  Following 
development,  YBa2Cu307_^  samples  were  etched  in  Nal/I+  filtered  semiconductor 
grade  gold  etch  for  45  s  and  rinsed  in  isopropyl  alcohol.  Cryogenic  ion  milling  for 
5  minutes  was  used  for  MgB2  samples  because  the  gold  etch  degraded  the  films. 

Fabricating  contact  layers  using  this  technique  was  superior  to  the  lift¬ 
off  contacts  for  many  reasons.  Having  a  noble  metal  covering  the  superconductor 
in  the  early  processing  steps  prevents  degradation  from  contact  with  water  and 
process  chemicals.  It  also  protects  from  any  thermal  damage  inflicted  by  the  ion 
mill.  It  lowers  the  contact  resistance  over  4  orders  of  magnitude  and  substantially 
improves  contact  adhesion 


■  YBCO 
□  Gold/YBCO 

IH  Alignment  Marks  (Gold/YBCO) 

Figure  3.5:  Central  region  of  the  sixteen  bridge  photolithographic  graph  mask  dark  lines  showing  the 
region  were  the  contact  was  removed  via  chemical  or  dry  etching. 


3.6  Tri-layer  Structure 

The  junctions  were  defined  in  the  bridges  by  using  electron  beam  litho¬ 
graphy  and  ion  damage.  Electron  beam  resist  is  not  thick  enough  to  serve  as  an 
implantation  mask  so  the  EBL  pattern  was  transferred  into  a  layer  of  germanium 
and  photoresist.  Samples  were  cleaned  in  isopropyl  and  dried  in  nitrogen  prior  to 
the  creation  of  a  tri-layer  structure.  (Figure  3.6).  They  were  spin-coated  with  an 


Figure  3.6:  Cross-section  showing  the  components  of  the  tri-layer  implant  mask. 
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800  nm  thick  layer  of  Shipley  S1808  photoresist  spun  at  4000  rpm  for  60  s,  fol¬ 
lowed  by  a  hard-bake  on  a  hot  plate  for  10  minutes  at  130°C.  This  was  followed  by 
thermal  evaporation  of  a  25  nm  thick  layer  of  amorphous  germanium.  A  100  nm 
layer  of  electron  beam  resist,  2%  polymethylmethacrylate  (PMMA)  in  chloroben¬ 
zene  was  spin-coated  on  top  of  the  Ge  at  4000  rpm  for  60  s.  The  sample  was  then 
baked  again  on  the  hot  plate  at  125°C  for  10  minutes.  It  was  important  for  the 
temperature  of  the  second  bake  to  be  lower  than  that  of  the  first,  otherwise  the 
underlying  photoresist  would  contract  and  crack  the  germanium  layer  as  shown  in 
Figure  3.7. 


Figure  3.7:  Microscope  picture  showing  a  cracked  tri-layer  mask.  The  circuit  pattern  is  visible  underneath 
the  cracked  resist/Ge  layers. 


3.7  Electron  Beam  Lithography 

For  positive  electron  beam  lithography  a  resist,  such  as  PMMA,  consisting 
of  long  molecular  chains  is  exposed  to  an  electron  beam  that  breaks  many  bonds 
in  the  chains.  Shorter  chains  are  more  soluble  in  developer  than  unbroken  chains. 
The  dose  of  electrons  required  to  break  enough  chains  to  make  the  exposed  area 
more  soluble  is  called  the  threshold  level  of  the  resist.  Thus  by  optimizing  the 
resist  with  the  proper  dose  of  electrons  and  developer,  fine  lithography  can  be 
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performed.  Features  of  the  order  of  10  nm  can  be  readily  obtained  in  PMMA 
using  commercial  high  voltage  lithography  systems. 

For  this  work  electron  beam  lithography  was  performed  on  a  JEOL  JSM- 
6400  scanning  electron  microscope  equipped  with  Nanometer  Pattern  Generation 
System  (NPGS)  to  modify  it  for  lithography  [37].  It  was  also  equipped  with  an 
electrostatic  beam  blanker  which  consisted  of  two  plates  in  the  column  that  the 
beam  passed  through.  A  potential  applied  across  the  plates  deflects  the  beam  so 
it  misses  the  aperture  and  prevents  the  beam  form  reaching  the  sample.  This  can 
be  used  to  prevent  exposure  of  critical  areas  while  performing  lithography. 

With  the  beam  optimized  the  sample  was  ready  to  be  positioned  for 
writing.  One  of  the  4  alignment  bars  next  to  the  common  electrodes  ( V~ ,  I~)  of  the 
photolithographic  pattern  (Figure  3.3)  was  located.  Magnification  was  increased 
so  that  the  bar  filled  the  whole  field  of  view  (10,000X).  This  limited  exposure  to 
the  resist  directly  above  the  bar  and  prevented  exposure  where  the  pattern  was  to 
be  written.  The  sample  was  then  moved  into  position  by  following  the  bar  into 
the  center.  At  the  end  of  the  alignment  bar  another  hole  was  burnt  in  the  PMMA 
to  adjust  focus. 

The  beam  was  then  blanked,  and  the  scan  coils  were  placed  under  the 
control  of  the  NPGS.  Magnification  was  changed  to  1000X  which  has  a  90  fj,m  x  90 
/jm  field  of  view.  A  program  for  alignment  was  then  executed.  The  program  scans 
4  rectangular  regions  on  the  outside  of  the  pattern  without  exposing  the  center. 
The  rectangles  are  then  moved  to  correspond  to  the  ends  of  the  4  alignment  bars  on 
the  sample.  The  program  then  calculates  a  transformation  matrix  to  correct  any 
x  and  y  error.  After  alignment  is  complete  the  run  file  is  called  and  the  DesignCad 
pattern  is  written.  Typical  run  files  were  sets  of  lines  used  to  define  the  junction 
barriers.  A  single  line  would  cross  8  of  the  16  bridges.  Figure  3.8  shows  the  pattern 
written  for  5-junction  arrays. 

Exposed  samples  were  then  developed  for  35  s  in  1:3  Methyl  Isobutyl 
Ketone  in  Isopropyl  Alcohol  (MIBK/IPA).  During  development  the  sample  was 


Figure  3.8:  Central  region  of  the  16-Bridge  mask  showing  5  lines  (junctions)  scribed  across  each  bridge. 


agitated  gently  to  promote  removal  of  the  exposed  region.  After  development 
samples  were  rinsed  in  isopropyl  and  blown  dried  with  compressed  nitrogen. 

Using  a  20  pA  beam  current  line  doses  between  3-8  nC/cm  were  used  to 
create  line  widths  ranging  between  30-100  nm.  The  dose  required  to  generate  a 
desired  line  width  varied  between  samples  and  had  to  be  optimized. 

3.8  Reactive  Ion  Etching 

Reactive  ion  etching  is  a  common  etching  technique  used  in  microelec¬ 
tronics  processing.  It  combines  physical  sputtering  and  chemical  etching  and  can 
be  used  to  anisotropically  etch  many  materials.  A  plasma  discharge  is  created 
in  the  presence  of  a  low  pressure  gas.  The  species  in  the  gas  disassociate  in  the 
plasma  and  form  chemically  reactive  radicals.  They  are  accelerated  to  the  tar¬ 
get  where  they  interact  chemically  to  form  gas  molecules  with  the  target  atoms 
that  are  then  pumped  away.  The  acceleration  makes  the  etch  highly  directional. 
Choosing  a  gas  whose  radicals  react  with  the  target  but  not  the  mask  results  in  a 
highly  anisotropic  etch  capable  of  producing  nearly  vertical  sidewalls. 

Developed  samples  were  loaded  into  a  Plasma  Science  table  top,  parallel 
plate  reactive  ion  etcher  (RIE).  This  system  had  a  turbo  pumped  8”  diameter 
cylindrical  quartz  chamber.  Inside  were  a  set  of  6”  circular  plates  with  2|”  spacing. 
The  lower  plate  was  cooled  with  chilled  water.  The  chamber  was  equipped  with 
two  process  gas  lines,  each  terminating  with  a  needle  valve  inline  with  a  shut  off 
valve  to  control  the  gas  flow.  One  line  could  be  used  to  introduce  either  Argon 
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or  Oxygen  and  the  other  was  used  for  CCI2F2,  CF3H,  CCI3F,  SFg  or  CF4.  A  300 
W  13.56  MHz  RF  power  was  used  to  generate  the  plasma  and  a  manual  matching 
network  was  used  to  match  the  impedance  of  the  plasma  to  the  power  supply. 

The  pattern  from  the  PMMA  was  transferred  into  the  Ge  using  a  CCI2F2 
plasma.  The  gas  pressure  was  maintained  at  14  mTorr,  and  the  dc  self  bias  at  200 
V.  The  forward  power  resulting  was  20  W.  The  matching  network  was  adjusted  so 
that  the  power  reflected  was  0  W.  This  etched  the  Ge  at  about  15  nm/min  and 
the  PMMA  at  about  30  nm/min.  Etching  was  performed  for  2  min.  to  ensure  the 
25  nm  layer  was  completely  etched.  Although  the  selectivity  between  PMMA  and 
Ge  was  not  very  good  for  this  etch  the  thickness  of  the  PMMA  was  sufficient  to 
transfer  the  pattern  without  broadening  of  the  features.  The  chamber  was  pumped 
back  out  to  its  base  pressure  and  the  sample  was  then  etched  in  Oxygen  plasma. 
The  gas  pressure  was  maintained  at  5  mTorr,  and  the  dc  self  bias  at  300  V.  The 
forward  power  resulting  was  28  W.  This  etched  the  photoresist  at  50  nm/min.  The 
sample  was  over  etched  for  1  min  to  ensure  that  the  layer  was  completely  etched. 
This  etch  also  removed  the  PMMA  layer  but  had  no  effect  on  Ge  or  YBCO. 

3.9  Ion  Damage 

To  create  the  actual  junction  the  masked  bridges  were  bombarded  with 
ions  in  an  ion  implanter.  The  ion  energy  was  balanced  so  that  it  was  high  enough 
to  penetrate  the  whole  thickness  of  the  unmasked  YBCO  film  to  create  uniform 
damage,  but  not  high  enough  to  penetrate  the  Ge/Resist  mask  and  damage  the 
electrodes.  Building  on  previous  work  the  majority  of  the  samples  were  implanted 
with  200  keV  Neon  with  doses  ranging  between  (5  x  1012  -  3  x  1013  ions/cm2). 

Etched  samples  were  attached  to  a  4”  wafer  with  silver  paint  and  sent  to  a 
company  for  ion  implantation  (Core  Systems)  [38].  To  limit  charging  the  surface 
of  the  sample,  it  was  grounded  to  the  wafer  with  silver  paint.  To  prevent  ion 
channeling  along  the  c-axis  of  the  crystal,  samples  were  tilted  at  7°  (Figure  3.10). 


Figure  3.9:  Grass  after  RIE  etching  of  the  tri-layer. 


Channeling  limits  the  amount  of  damage  because  the  ions  can  travel  very  far 
through  the  sample  before  experiencing  a  collision.  Because  of  the  high  aspect 
ratio  of  our  mask  it  was  important  to  ensure  that  the  tilt  angle  was  along  the  slits 
used  to  define  the  junctions.  Any  misalignment  would  create  a  shadow  and  result 
in  reduced  or  no  damage  depending  on  the  level  of  misorientation.  A  4”  copper 
wafer  was  machined  with  a  7°  inclined  ramp  in  the  center  (Figure  3.11).  This 
allows  us  to  align  the  sample  with  regard  to  the  tilt  rather  than  relying  on  the 
technicians  at  the  company.  The  company  could  then  implant  with  any  rotation 
at  a  0°  wafer  tilt  and  the  lines  were  always  7°  to  the  ion  beam.  The  machine  used 
was  a  Varian  CF  3030  implanter.  The  ion  beam  was  electrostatically  scanned. 
The  horizontal  and  vertical  scan  rates  were  256  Hz  and  2.250  KHz  respectively. 
The  machine  scans  an  area  of  a  4”  wafer  plus  a  of  over  scan  to  ensure  uniform 
implantation.  Uniformity  studies  by  Core  show  the  uniformity  of  the  implant  to  be 
between  0. 5-1.0%.  Table  3.1  shows  the  different  implantations  performed.  Most  of 
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Figure  3.10:  Etched  tri- layer  mask  depicting  the  7°  tilt  used  to  prevent  channeling. 

the  implantations  were  200  keV  Ne,  however  Ag  was  also  attempted.  The  current 
is  the  beam  current  used  over  an  area  of  a  4”  wafer.  Holding  dose  constant  a 
decrease  in  current  increases  implantation  time.  This  was  varied  to  investigate  the 
possibility  of  charging  effects. 


Figure  3.11:  Sample  holder  designed  to  prevent  misalignment.  The  center  of  the  wafer  contains  a  ramp 
with  a  7°  tilt.  The  lines  scribed  by  EBL  are  alighned  parallel  with  the  ramp  and  the  wafer  is  implanted 
at  0°. 
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Table  3.1:  Implantations 


Energy(keV) 

Ion 

Dose(^f) 

Current  (fJ-A) 

100 

Ne 

1.0  x  1013 

1.2 

200 

Ne 

5.0  x  1012 

1.0 

200 

Ne 

1.0  x  1013 

0.12 

200 

Ne 

1.0  x  1013 

0.5 

200 

Ne 

1.0  x  1013 

1.0 

200 

Ne 

1.0  x  1013 

1.2 

200 

Ne 

1.5  x  1013 

0.5 

200 

Ne 

1.5  x  1013 

1.0 

200 

Ne 

2.0  x  1013 

0.5 

200 

Ne 

2.0  x  1013 

1.0 

200 

Ne 

3.0  x  1013 

1.0 

400 

Ag 

1.0  x  1011 

0.32 

400 

Ag 

1.0  x  1012 

1.0 

400 

Ag 

1.0  x  1013 

1.0 

3.9.1  Simulations 

Simulations  were  performed  using  SRIM  2003  (Stopping  and  Range  of 
Ions  in  Matter)  to  help  optimize  the  implantation  process  [39].  SRIM  is  a  freeware 
program  used  in  ion  implantation  to  calculate  the  stopping  and  range  of  ions 
in  matter  using  a  quantum  mechanical  treatment  of  ion-atom  collisions.  SRIM 
simulations  showed  that  for  a  200  keV  Ne+  implantation  the  ions  are  eneregetic 
enough  to  penetrate  an  unmasked  200  nm  thick  YBCO  film  (junction  barriers)  and 
that  the  Ge/photoresist  mask  was  of  sufficient  density  to  stop  ions  from  damaging 
the  electrodes  (areas  underneath  the  trilayer  mask). 

SRIM  is  a  useful  tool  for  investigating  the  range  in  different  targets,  but 
it  is  difficult  to  do  2-dimensional  device  modeling.  It  also  lacks  the  capability 
to  take  into  account  the  crystallinity  of  the  sample  where  channeling  effects  can 
significantly  alter  the  damage.  To  perform  more  complex  simulations  a  commercial 
software  package  called  Athena  from  Silvaco  [40]  was  used.  It  is  a  common  tool 
used  in  semiconductor  processing  for  the  simulation  of  implants,  annealing  and 
etching.  The  model  used  for  the  simulations  was  a  Monte  Carlo  model  based  on 
the  Binary  Collision  Approximation. 
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Figure  3.12  is  an  example  of  a  simulation  that  shows  the  structure  and 
implantation  damage  from  lxlO13  Ne+/cm2.  The  mask  opening  is  parallel  to  the 
b-c  plane  of  YBa2Cu307_i  and  the  cross  sectional  view  is  taken  along  the  a—c  plane. 
Implantation  parameters  rotation,  tilt,  energy  and  dose  were  varied  along  with  the 
mask  dimensions  and  materials  to  gain  insight  into  this  process.  Simulations  were 
also  performed  for  comparison  with  experimental  results. 

3.10  Electron  Beam  Proximity  Effect 

In  EBL  backscattered  electrons  can  travel  tens  of  microns  in  the  resist. 
The  range  of  the  electron  beam  scattering  depends  on  the  energy  of  the  incident 
electron  beam,  the  thickness  of  the  resist,  and  the  materials  under  the  resist.  If 


Figure  3.12:  Implantation  simulation  of  a  400  nm  space  junction  pair.  The  contour  scale  is  the  base  ten 
log  of  the  number  of  displacements  for  cm3. 
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enough  of  these  electrons  expose  an  area  over  the  threshold  of  the  resist  it  can  result 
in  the  development  of  undesired  features.  This  effect  called  the  proximity  effect 
and  typically  manifests  itself  in  the  broadening  of  the  pattern  features.  Since  the 
density  of  the  scattered  electrons  decay  exponentially  with  distance  from  the  beam, 
closely  distributed  patterns  are  affected  by  the  proximity  effect  substantially  more 
than  sparsely  distributed  ones.  In  the  case  of  a  set  of  uniformly  spaced  identical 
parallel  lines  for  the  ion-damaged  Josephson  junction  arrays,  the  resultant  pattern 
will  not  faithfully  replicate  that  which  was  written.  Lines  in  the  center  of  the 
pattern  receive  a  higher  dose  than  those  on  the  outside  resulting  in  thinner  line 
widths  for  those  written  on  the  edges. 

Figure  3.13  shows  an  exposure  vs.  distance  curve  determined  by  measur¬ 
ing  the  line  widths  for  different  doses.  A  double-Gaussian  function  was  fit  to  the 
data  for  comparison.  A  double-Gaussian  model  is  commonly  accepted  as  an  ap¬ 
proximation  to  the  real  scattering  behavior.  This  curve  illustrates  how  important 
the  proximity  effect  is.  The  equivalent  relative  exposure  values  for  distances  of  40, 
250  and  500  nm  are  approximately  7  x  10~2,  4  x  10~3,  and  2  x  10~3,  respectively. 
This  means  that  writing  an  80  nm  wide  single  line  will  expose  the  area  250  nm  and 
500  nm  away  from  the  line  to  ^  and  ^  of  the  threshold  of  the  resist  required  for 
development.  This  may  not  seem  to  be  a  serious  problem  for  just  a  few  lines,  how¬ 
ever  as  the  number  of  lines  increases  or  the  distance  between  the  lines  decreases, 
the  exposure  contributions  from  scattering  of  all  the  other  line  writings  may  add 
up  to  a  significant  value  great  enough  to  impact  the  line  width.  If  the  same  dose 
is  used  for  writing  each  line,  the  line  at  the  middle  of  the  array  sees  more  exposure 
than  the  lines  at  the  far  ends. 

Attempts  to  modulate  the  dose  for  each  line  to  compensate  for  the  dif¬ 
ferences  a  line  in  the  middle  would  “see”  compared  with  one  at  the  edge  was 
complicated  and  did  not  provide  satisfactory  results.  The  compensation  was  ei¬ 
ther  too  much  or  not  enough.  We  attributed  this  failure  to  the  inadequacy  of  the 
double-Gaussian  approximation,  which  may  be  too  simple  to  properly  take  into 
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Figure  3.13:  Exposure  vs.  distance  data  for  single  line  writing  demonstrating  the  electron  beam  proximity 
effect.  Distance  is  measured  from  the  center  of  the  line.  The  solid  line  is  a  double  Gaussian  fit. 

account  the  real  behavior  of  electron  scattering  in  our  samples. 

Instead  of  modulating  the  line  doses,  the  dose  was  offset  by  writing  ad¬ 
ditional  simple  patterns  outside  the  line  array.  Writing  additional  lines  with  the 
same  dose  would  make  the  lines  in  the  middle  part  of  the  array  experience  a  more 
uniform  amount  of  electron  scattering  from  the  surroundings.  However,  the  addi¬ 
tional  lines  would  then  develop  and  be  part  of  the  pattern.  Instead,  the  additional 
lines  can  be  replaced  with  simple  rectangles  with  the  same  average  surface  dose. 
Since  most  of  the  areas  of  the  rectangles  are  far  away  from  all  the  lines,  the  lines  in 
the  array  still  feel  roughly  the  same  amount  of  scattering  from  writing  the  rectan¬ 
gles  as  if  they  were  lines.  The  writing  of  the  rectangles  did  not  leave  any  pattern 
after  developing  as  long  as  the  dose  was  below  the  threshold  of  the  resist.  Ex¬ 
periments  have  shown  that  this  idea  works  very  well.  Figure  3.14  illustrates  the 
writing  pattern  for  100  parallel  lines.  Each  line  is  80  [im  long  and  was  written 
by  scanning  the  electron  beam  in  only  one  scan.  The  width  is  controlled  by  the 
writing  dosage.  The  center-to-center  distance  between  adjacent  lines  was  250  nm 
and  500  nm  for  two  separate  samples.  Two  rectangles  were  written  next  to  the 
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line  at  the  each  end  of  the  array  as  compensation  for  the  proximity  effect.  Their 
size  was  80  pm  by  20  pm.  The  separation  between  the  rectangle  and  the  last  line 
in  the  array  was  kept  the  same  as  the  line  period  in  the  array.  The  doses  for  all  the 
lines  in  an  array  are  the  same  and  the  doses  for  the  two  corresponding  rectangles 
are  the  same.  The  doses  for  the  rectangles  are  set  by  the  average  surface  dose  of 
the  arrays.  In  order  to  determine  the  effect  of  the  rectangles,  another  pattern  with 
a  rectangle  at  only  one  end  was  also  written.  The  doses  for  writing  the  rectangles 
are  80  and  120  pC/cm2  for  250  and  500  nm  period  line  arrays.  The  dose  threshold 
for  developing  PMMA  is  about  250  pC/cm2.  Consequently,  writing  compensation 
rectangles  left  no  pattern  after  development.  After  developing  and  RIE  etching, 
the  stencil  structure  was  measured  under  the  same  SEM.  The  NPGS  was  used  to 
record  the  SEM  images  and  measure  the  widths  of  the  resultant  lines. 

Figure  3.15  (a)(b)(c)  shows  the  left  end,  middle  and  the  right  end  region 
of  the  100-line  (spaced  at  500  nm)  array  with  compensation  rectangles;  (d)(e)(f) 
shows  the  corresponding  regions  with  a  compensation  rectangle  written  only  at 
the  right  end.  The  widths  of  the  lines  for  the  first  sample  are  all  86.9  nm,  while 
the  widths  of  the  lines  for  the  second  one  are  80.7,  93.1  and  93.1  nm  for  the  left, 


100  lines 


Compensation  Rectangles 


Figure  3.14:  Illustration  of  the  pattern  for  writing  100  parallel  lines  and  compensation  rectangles  at  each 
end  of  the  lines  (not  to  scale). 
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(d)  (e)  (0 

Figure  3.15:  SEM  pictures  of  the  100-line  array  with  500  nm  period,  (a),  (b),  (c)  The  pictures  of  the 
left  end,  middle,  and  the  right  end  of  the  line  array  of  a  sample  with  both  compensation  rectangles 
written,  (d),  (e),  (f)  The  pictures  of  the  corresponding  parts  of  the  line  array  of  a  sample  with  only  right 
compensation  triangle  written. 

middle,  and  right  region.  It  is  clear  that  the  line  width  differences  induced  by  the 
electron  proximity  effect  have  been  corrected  with  the  rectangles.  The  result  for  a 
sample  of  100-line  with  250  nm  period  is  shown  in  Figure  3.16.  The  array  has  the 
same  line  width  37.2  nm  for  lines  in  the  middle  and  by  the  edge  of  the  array. 

The  sizes  of  the  compensation  rectangles  are  not  critical.  The  compen¬ 
sation  rectangles  only  have  to  cover  the  entire  area  from  which  scattered  electrons 
could  reach  the  array.  From  Figure  3.13,  the  range  of  electron  scattering  is  less 
than  10  fim.  Therefore,  the  size  of  the  rectangles  were  chosen  to  be  80  /Ltm  by 
20  fim.  Since  rectangles  were  used,  which  are  continuous  patterns,  to  compensate 
scattering  from  discrete  lines  (or  individual  patterns  in  two  dimensional  cases), 
summation  is  substituted  by  an  integral  in  the  approximation,  which  is  good  only 
when  the  array  size  is  sufficiently  large  and  the  period  is  small  enough.  Based  on 
Gaussian-like  behavior  of  electron  scattering,  it  is  estimated  that  the  array  size  has 
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(a) 


(b) 


Figure  3.16:  SEM  pictures  of  the  100-line  array  with  250  nm  period,  (a)  The  middle  part,  (b)  the  right 
end  of  the  line  array.  A  compensation  rectangle  has  been  written  by  the  right  end. 


to  be  larger  than  the  electron  scattering  range.  From  experiments,  this  method 
works  with  periods  of  500  nm  and  less. 


Chapter  4 


YBa2Cii307_£  Josepshson 
Junctions 

4.1  Single  Junctions 

A  successful  VLSI  junction  technology  must  be  able  to  fabricate  junc¬ 
tions  with  junction  parameters  Ic  and  Rn  that  do  not  to  vary  more  than  10% 
between  different  junctions.  Furthermore,  different  applications  require  different 
values  of  IcRn.  For  the  Josephson  arbitrary  waveform  generator,  an  IcRn  product 
between  24  and  40  fiV  at  65  K  is  desired,  however  SQUID  applications  require  an 
IcRn  >  100  jiV .  Therefore  not  only  must  a  successful  technology  create  uniform 
parameters,  it  must  also  be  able  to  reproduce  a  given  value  of  IcRn. 

Although  the  primary  goal  of  this  work  was  to  fabricate  junction  arrays, 
single  junctions  were  also  studied  so  that  individual  junction  properties  could  be 
compared  with  those  of  arrays.  Furthermore  the  uniformity  of  junction  parameters 
is  difficult  to  determine  in  an  array  because  it  would  require  voltage  probes  in 
between  each  of  the  junctions  which  is  not  practical. 

Current  voltage  {I —  V)  measurements  were  performed  to  determine  the 
uniformity  and  reproducibility  of  the  junction  parameters  Ic  and  Rn.  I  —  V  char¬ 
acteristics  could  be  classified  into  two  types.  Restively  shunted  junctions  similar 
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to  those  described  by  the  Stewart-McCumber  model  [11,  12]  (RSJ)  and  curves 
rounded  near  Ic  resembling  those  of  flux-flow  “long”  junctions. 

Figure  4.1  shows  the  I—V  characteristic  at  63  K,  for  a  junction  fabricated 
from  200  nm  YE^CusOr-a  on  LaAlOs,  implanted  with  2xl013  Ne+/cm2  at  200 
keV.  Figure  4.1(a)  is  a  fit  to  the  data  using  the  RSJ  model  (equation  2.17).  The 
data  certainly  follows  the  trend  of  the  RSJ  model  however,  it  does  not  fit  well  near 
Ic  where  the  critical  current  is  rounded.  At  higher  biases  the  resistance  doesn’t 
follow  that  of  the  RSJ  model  due  to  broadening  of  the  barrier.  Rounding  of  the 
I  —  V  characteristic  near  the  critical  current  can  occur  in  junctions  from  thermal 
fluctuations  that  disrupt  the  coupling  of  the  phases  of  the  order  parameters  of 
the  two  superconducting  electrodes.  A  common  method  for  modeling  of  thermal 
rounding  in  Josepshon  junctions  was  first  described  by  Ambegaokar  and  Halerpin 
[41].  A  thermal  noise  current  was  introduced  as  an  additional  term  in  the  RSJ 
model  to  account  for  the  thermally  activated  phase  slip  processes.  This  model  was 
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Figure  4.1:  I  —  V  characteristic  at  63  K  showing  a  fit  with  (a)the  RSJ  model,  and  (b)  the  RSJ  model 
with  thermal  rounding  and  a  current  dependant  resistance. 
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combined  with  the  current  dependant  resistance  used  by  Katz  and  Dynes  (equation 
2.35)  to  fit  the  data.  The  result  is  shown  in  Figure  4.1(b).  The  noise  temperature 
calculated  from  the  fit  was  71.9  K,  which  is  close  to  the  measurement  temperature 
of  63  K.  The  additional  noise  temperature  is  attributed  to  the  measurement  system. 
Both  of  the  fits  produce  similar  IcRn,  35  fiV,  and  31  fiV  for  the  RSJ  and  modified 
RSJ  models  respectively. 

Figure  4.2  shows  the  I  —  V  characteristics  at  5  different  temperatures, 
for  a  junction  fabricated  from  200  nm  YBa2Cu307_a  on  LSAT,  implanted  with 
lxlO13  Ne+/cm2  at  200  keV.  Near  the  Tc  of  the  junction,  I  —  V  characteristics 
have  a  RSJ-like  shape.  However  as  temperature  is  reduced  the  junction  becomes 
strongly  coupled  and  the  I  —  V  curve  becomes  rounded  near  Ic.  These  rounded 
curves  are  similar  to  the  rounding  seen  in  “long”  junctions,  where  the  bridge  width, 
w  is  larger  then  the  Josephson  penetration  depth,  w/Xj.  The  crossover  from  short 
to  long  junction  behavior  is  generally  taken  to  be  when  the  ratio  w/Xj  >  6.  This 
was  calculated  for  the  temperatures  of  the  I—V  curves  in  Figure  4.2  and  is  displayed 
in  the  legend.  At  the  crossover  temperature,  T  ~  76.9  K,  wjXj  is  over  a  factor  of 
3  smaller  than  what  would  be  expected  if  this  was  solely  a  long  junction  effect. 


Voltage  (nV) 


Figure  4.2:  I  —  V  characteristics  measured  for  the  same  junction  at  five  temperatures.  From  top  to 
bottom  :  75.6  K,  76.9  K,  77.8  K,  78.5  K,  79.7  K.  The  junction’s  barrier  undergoes  a  transition  from  weak 
coupling  to  strong  coupling  with  decreasing  temperature. 
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Figure  4.3  shows  the  flux-flow  regime  I  —  V  characteristics  of  the  ion 
damage  junction  from  Figure  4.2.  They  have  V  oc  7n  with  n  increasing  with 
current  density.  This  is  consistent  with  models  of  vortex  motion  used  to  describe 
non-linear  I —  V  characteristics  in  flux-flow  junctions  [9,42].  The  transition  from 
RSJ-like  to  flux-flow  can  be  attributed  to  the  junction  changing  from  a  SN'S  to 
a  SS'S  junction  as  temperature  is  reduced  below  that  of  the  critical  temperature 
of  the  weaklink.  From  Likharev’s  theory  of  SS’S  junctions  [9]  the  current  phase 
relationship  becomes  multi-valued  when, 

m  z  3.5 

after  which  Abrikosov  vortices  may  enter  the  weaklink.  The  entrance  and  motion  of 
these  vortices  determine  the  properties  of  the  junction,  some  of  which  resemble  the 
ideal  Josephson  effects.  When  T  — >  Tc',  the  critical  temperature  of  the  weaklink, 
f  — *  oo  and  L/£  — >  0.  Thus  ideal  Josephson  effects  are  observed  near  T'c  for  large 
barriers.  Estimates  of  the  coherence  length  using  the  Likharev  model  at  the  cross¬ 
over  temperature  are  consistent  with  values  determined  from  proximity  effect  fits 
of  the  critical  current  temperature  dependance  [43]. 


Figure  4.3:  Flux-flow  V— I  characteristics  from  the  same  junction  in  Figure  (4.2).  The  solid  lines  are  fits 
of  the  form  V  oc  In,  where  the  exponent  n  is  related  to  vortex  motion  in  the  weaklink. 
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Figure  4.4:  Critical  current  and  normal  resistance  vs.  temperature  for  the  same  junction  from  Figure 
4.2.  The  inset  shows  the  IcRn  product  as  a  function  of  temperature.  It  has  a  peak  where  the  junction 
transitions  from  RSJ-like  to  strongly  coupled. 


Figure  4.4  shows  the  temperature  dependance  of  the  critical  current  and 
normal  state  resistance  for  the  same  junction  in  Figure  4.2.  The  inset  shows 
that  IcRn  has  a  maximum  near  the  crossover  temperature.  This  feature  has  been 
empirically  observed  in  many  junctions  near  the  temperature  where  the  junction 
becomes  strongly  coupled. 

The  uniformity  of  the  junction  parameters  for  the  single  junctions  fabri¬ 
cated  on  this  sample  were  outstanding.  Figure  4.5  shows  the  I—V  characteristics 
measured  at  79.7  K,  for  6  of  8  single  junctions  fabricated  on  one  side  of  a  16-bridge 
sample  by  writing  a  single  line  across  all  bridges.  Two  of  the  junctions  were  not 
measured  due  to  faulty  contacts.  The  junctions  are  labeled  J1-J8  with  the  number 
representing  the  position  of  the  bridge.  The  I  —  V  curves  for  J3,  J4,  J6,  and  J7 
are  nearly  identical,  which  seldom  occurs  in  other  HTS  technologies.  J1  and  J8 
have  larger  critical  currents  and  smaller  resistances  indicating  less  damage.  These 
junctions  are  located  on  the  top  and  bottom  bridges  on  the  outside  of  the  EBL 
pattern.  Junctions  fabricated  on  these  bridges  consistently  had  less  damage  for 
many  samples.  It  was  suspected  that  the  substrate  on  the  outside  of  the  bridges 
was  charging  during  the  implantation  and  repelling  incoming  ions  and  hence  reduc¬ 
ing  the  damage.  The  implantation  time  was  increased  by  an  order  of  magnitude 
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Figure  4.5:  Six  junctions  from  the  same  chip.  The  legend  labels  the  curves  from  top  to  bottom.  The 
parameters  of  J7,  J3,  J4,  and  J6  are  nearly  identical. 

by  decreasing  the  current  to  allow  time  for  accumulated  charge  to  dissipate.  This 
had  no  effect  and  the  outside  bridges  still  appeared  less  damaged  then  those  in  the 
center. 

The  temperature  dependance  of  the  junction’s  critical  current  was  con¬ 
sistent  with  the  Katz-Dynes  model  which  suggests  a  (1  -  T/Tc )3  dependance  near 
Tc  and  a  (1  -  T/Tc )2  dependance  at  lower  temperatures.  Figure  4.6  shows  data 
for  a  a  junction  fabricated  from  100  nm  YBa2Cu307_^  on  AI2O3,  implanted  with 
lxlO13  Ne+/cm2  at  200  keV,  showing  both  fits. 
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Figure  4.6:  Critical  current  vs.  temperature,  near  Tc  the  data  have  a  (1  —  T/Tc)3  dependance.  At  lower 
temperatures  it  follows  a  (1  —  T/Tc)2  dependance. 
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Figure  4.7:  Shapiro  steps  in  a  junction’s  I—V  characteristics  for  different  applied  10  GHz  RF  powers, 
measured  at  77.7  K.  The  y-axis  has  been  offset  for  clarity. 


RSJ-like  junctions  exhibited  the  ac  Josepshon  effect  [8].  Steps  occurred 
precisely  at  V  —  with  harmonics  out  to  n  >  10.  Figure  4.7  shows  I  — 
V  characteristics  for  J6  of  Figure  4.5  measured  at  77.7  K.  10  GHz  microwave 
radiation  was  introduced  through  an  antenna  and  curves  were  measured  for  6 
different  power  levels.  The  amplitude  dependance  of  the  Shapiro  steps  as  a  function 
of  RF  Voltage  for  another  similar  junction  was  measured  and  is  shown  in  Figure 
4.8.  The  amplitudes  have  a  Bessel-like  dependance  confirming  that  the  steps  are 
due  to  the  ac  Josephson  effect. 


Figure  4.8:  Bessel-like  dependencies  of  the  critical  current  (n  =  0),  and  the  first  two  Shapiro  steps.  The 
inset  shows  the  full  range  of  the  data. 
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The  modulation  of  the  supercurrent  in  magnetic  field  was  also  measured 
as  a  test  for  the  dc  Josephson  effect.  Figure  4.9  shows  a  junction  fabricated  from 
200  nm  YBa2Cu307_5  on  LaAlC>3,  implanted  with  1.5xl013  Ne+/cm2  at  200  keV. 
The  slit  in  the  mask  used  to  define  the  ion  damage  region  was  50  nm.  The  data  were 
taken  at  65  I<  by  sweeping  the  field  from  -40  to  40  Gauss  and  then  sweeping  back. 
The  two  traces  represent  data  taken  with  increasing  and  decreasing  field.  Hystere¬ 
sis  is  apparent  between  the  two  sweeps.  This  occurs  when  the  field  exceeds  Hc i 
and  Abrikosov  vortices  enter  the  film  and  become  trapped  in  the  electrodes  near 
the  junction.  The  vortices  partially  cancel  the  Meissner  screening  currents  near 
the  edge  of  the  junction  and  shifts  the  pattern  [44],  The  period  of  magnetic  field 
oscillation  was  calculated  to  be  15.6  Gauss  by  using  equation(2.26),  and  estimating 
A l  —  300  nm.  This  is  in  close  agreement  to  the  measured  value  of  the  average 
oscillation  period,  14.9  Gauss.  This  is  a  rare  occurrence  in  Josephson  junctions 
with  a  planar  geometry.  Planar  junctions  usually  exhibit  much  smaller  oscillation 
periods  attributed  to  flux  focusing  effects.  Theories  by  Rosenthal  et.  al.  [45],  and 
Humphries-Edwards  [46]  are  commonly  used  to  predict  the  periods  in  planar  junc¬ 
tions.  They  make  the  claim  that  in  a  planar  geometry  the  junction  is  unshielded 
opposed  to  the  vertical  geometry  in  classical  junctions  with  thick  electrodes.  More 


Figure  4.9:  Fraunhofer  pattern  demonstrating  the  dc  Josephson  effect.  The  two  curves  show  the  hysteresis 
for  sweeps  of  increasing  and  decreasing  field. 
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field  can  penetrate  in  the  planar  junction  case  giving  the  appearance  of  a  larger 
junction  area.  The  period  of  oscillation  for  the  junction  in  Figure  4.9  was  cal¬ 
culated  using  these  theories,  resulting  in  9.5  and  4.7  Gauss  for  Rosenthal,  and 
Humpries-Edwards  respectively.  Both  of  these  theories  accurately  predict  the  os¬ 
cillation  periods  in  grain  boundary  and  ramp  junctions  yet  they  clearly  fall  short 
in  the  case  of  the  ion  damage  junctions  where  the  classical  model  describes  the 
data  better. 

Figure  4.10  shows  the  I—V  characteristic  without  field  and  in  field  at  the 
first  minimum  (<F0).  A  non-Josephson  current  of  25  fiA  is  observed  and  a  common 
occurrence  in  ion-damage  YBa2Cu307_^  junctions  [18,47].  It  can  occur  from  a 
number  of  different  mechanisms  that  are  related  to  strong  coupling  effects  [48,49]. 
It  is  insensitive  to  the  phase  difference  across  the  junction  and  does  not  modulate 
in  magnetic  field.  It  can  be  viewed  as  a  superconducting  short  in  parallel  with  the 
junction.  This  is  consistent  with  I—V  temperature  data  (figure  (4.2)  and  simulation 
results  that  show  the  barrier  is  made  up  of  grains  of  different  damage  levels.  As 
temperature  decreases  some  of  the  grains  become  superconducting  and  create  a 
percolative  network  of  superconducting  shorts.  Non-uniform  current  distribution 
across  the  barrier  from  defects  can  also  prevent  the  minima  from  reaching  zero 


Figure  4.10:  I  —  V  characteristic  with  and  without  an  applied  flux  of  4>o. 
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thus  the  value  of  the  current  at  the  minima  can  only  be  taken  as  an  upper  limit 
on  the  excess  current  in  the  junction. 

The  effects  of  anneals  on  the  junctions  was  also  investigated  in  an  attempt 
to  improve  uniformity.  Anneals  were  performed  in  Oxygen  at  800,  and  400  °C  and 
also  in  air  at  100°C.  The  two  highest  temperature  anneals  resulted  in  reduced  Tc  of 
the  bulk  films  and  in  some  cases  insulating  samples.  The  low  temperature  anneal 
on  the  other  hand,  substantially  improved  the  devices. 

Figure  4.11  shows  the  IcRn  vs.  temperature  for  the  same  junction  before 
and  after  a  2  hour  100°C  anneal.  The  operating  temperature  of  the  junction 
increased  by  over  ten  degrees.  IcRn  also  showed  a  considerable  increase  due  to  an 
increase  in  critical  current.  Figure  4.12  shows  a  plot  of  IcRn  vs  Rn,  for  Rn  of  0.07 
and  0.05  there  is  a  5  times  and  2  times  increase  in  IcRn  respectively.  Figure  4.13 
compares  two  I  —  V  curves  of  the  same  resistance  showing  that  it  retains  its  RSJ- 
like  appearance  after  a  2  times  increase  in  the  critical  current.  This  is  significant 
because  high  IcRn  products  are  desirable  for  some  devices  such  as  SQUIDs  or 
SQIFs. 

It  seems  likely  that  at  100°C  the  thermal  energy  is  sufficient  for  the 
recombination  of  vacancy-interstitial  defects  which  results  in  recovery  of  Tc.  At 


Figure  4.11:  IcRn  vs.  temperature  for  the  same  junction  before  and  after  a  2  hour  anneal  at  100°C. 
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Figure  4.12:  IcRn  vs.  Rn  showing  a  substantial  increase  in  IcRn,  for  constant  Rn. 


the  higher  temperatures  the  thermal  energy  is  greater  then  the  oxygen  activation 
energy  which  results  in  oxygen  loss  and  a  lower  Tc. 

The  results  from  the  studies  of  the  single  junctions  demonstrate  that 
this  technology  is  capable  of  creating  uniform  and  reproducible  HTS  Josepshon 
junctions  with  RSJ-like  I  —  V  characteristics.  Furthermore  the  anneals  provide  a 
way  to  substantially  increase  IcRn  which  broadens  the  possible  applications. 


Figure  4.13:  I—V  characteristics  before  and  after  the  anneal.  The  resistance  is  the  same  and  the  critical 
current  is  increased  by  a  factor  of  two.  After  the  anneal  the  I  —  V  characteristic  retains  its  RSJ-like 
shape. 
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4.2  Junction  Pairs 

Most  applications  associated  with  Josephson  junction  arrays  require  that 
all  the  Josephson  junctions  in  the  array  are  phase-locked  to  the  applied  microwave 
or  to  each  other  [22].  The  interactions  between  Josephson  junctions  have  been 
investigated  intensively  [50] .  However,  due  to  the  complexity  of  Josephson  junc¬ 
tions  as  highly  non-linear  quantum  devices,  the  physics  of  the  interactions  has  not 
been  fully  understood.  Interaction  mechanisms  (e.g.,  mediated  by  Cooper  pairs, 
quasi-particles  and  electromagnetic  field)  have  been  suggested  to  account  for  the 
mutual  phase- locking  between  the  junctions  [50].  Because  of  the  short  coherence 
length  in  HTS  materials  the  interaction  mediated  by  Cooper  pairs  is  expected  to  be 
weak  unless  the  junctions  are  extremely  close  to  each  other.  Theoretical  analysis 
shows  that  apparent  mutual  interactions  can  be  observed  if  the  junction  distance 
is  less  than  50£  [51].  The  fabrication  difficulty  of  closely  spaced  junctions  has 
ruled  out  many  common  types  of  HTS  Josephson  junctions,  such  as  artificial  grain 
boundary  junctions  and  ramp- type  junctions.  Furthermore  it  is  a  difficult  task  to 
fabricate  pairs  with  an  electrode  in  between,  so  that  individual  junction  behavior 
can  be  measured.  EBL  was  used  to  pattern  two  slits  with  a  center  electrode  in 
the  tri-layer  mask  (Figure  4.14(a))  at  the  center  of  the  16  bridge  pattern  (Figure 
3.3)  where  two  of  the  bridges  intersect  the  common  electrode.  After  ion  damage 
the  uncovered  material  has  a  reduced  Tc  which  effectively  transfers  the  pattern  in 
to  the  film  with  no  need  for  etching.  Figure  4.14(b)  shows  the  measurement  con¬ 
figuration  for  simultaneous  measurement  of  individual  I —  V  characteristics  along 
with  the  series  pair. 

Many  identical  pairs  were  created  with  200,  400,  and  800  nm  center  to 
center  spacing.  Figure  4.15(a)  shows  I —  V  curves  measured  at  63  K  for  a  400 
nm  spaced  pair  with  slit  width  of  50  nm,  fabricated  from  a  200  nm  YBa2Cu307_<5 
film  on  LaC>3  implanted  with  2xl013  Ne+/cm2  at  200  keV.  The  analog  oscilloscope 
was  used  so  that  all  of  the  traces  could  be  viewed  simultaneously  because  the 
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Figure  4.14:  The  SEM  picture  of  a  junction  pair  after  ion  damage.  A,  B,  and  C  denote  the  photoresist 
mask,  the  YBa2Cu307_j  bridge,  and  the  sapphire  substrate,  respectively,  (b)  Schematic  drawing  of  the 
junction  pair  and  measurement  configuration. 


digital  oscilloscope  only  had  two  channels.  Traces  A  and  B  are  for  the  individual 
junctions  and  trace  C  is  for  the  pair  in  series.  Figure  4.15(b)  is  with  12  GHZ 
applied  microwave  radiation.  Large  amplitude  Shapiro  steps  are  visible  in  the 
single  junctions  at  V  =  —  n  x  24  \iV  and  twice  this  amount  in  the  series  pairs. 

It  is  assumed  that  these  junctions  are  locking  to  the  drive  signal  because  other 
experiments  did  not  show  any  evidence  of  an  order  parameter  interaction. 


(a)  (b) 


Figure  4.15:  Oscilloscope  pictures  of  the  current  voltage  characteristics  of  junction  1  (curve  A),  junction 
2  (curve  B,  and  the  pair  in  series  (curve  C)  without  (a)  and  with  (b)  12  GHz  microwave  radiation. 
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Figure  4.16  shows  results  from  a  150  nm  spaced  pair  fabricated  from  200 
nm  YBa2Cus07_5  on  LaA103  implanted  with  lxlO13  Ne+/cm2  at  200  keV.  The 
I  —  V  curves  A  and  B  are  measured  using  the  center  electrode,  and  across  the 
two  in  series,  C.  Under  13  GHz  the  series  pair  showed  steps  that  followed  the 
AC  Josephson  relationship  for  a  single  junction.  Fractional  Shapiro  steps  were 
observed  on  the  other  two  curves.  At  low  temperatures  these  junctions  show  two 
junction  behavior.  This  is  strong  evidence  that  the  damage  merged  underneath  the 
electrode  to  form  a  single  barrier  and  contradicts  a  report  that  the  barrier  length 
in  these  junctions  is  equivalent  to  the  slit  size  in  the  mask  [21],  This  geometry 
provides  a  novel  tool  to  probe  inside  the  barrier  of  a  Josephson  junction. 

Athena  ion  implantation  simulations  were  performed  for  the  device  of 
Figure  4.16.  The  structure  and  results  are  shown  in  Figure  4.17.  The  implant 
damage  straggles  underneath  the  center  mask  region  and  merges  into  one  another. 
The  central  region  appears  to  have  slightly  lower  damage  or  equivalently  a  higher 
Tc,  which  explains  why  they  demonstrate  double  junction  behavior  at  low  temper¬ 
atures. 

Nearly  identical,  closely  spaced  junction  pairs  were  fabricated  with  elec- 


Figure  4.16:  I —V  characteristics  at  71  K  for  a  150  nm  spaced  junction  pair  with  13  GHz  microwave 
radiation.  A  and  B  were  measured  using  the  central  electrode  and  C  is  measured  across  the  pair. 
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Figure  4.17:  Implantation  simulation 


trades  in  between  them  to  study  individual  and  collective  junction  behavior.  No  or¬ 
der  parameter  interactions  or  phase  locking  was  observed  possibly  due  to  the  short 
coherence  length  in  YE^CusOy-j.  Junctions  with  spacings  of  150  nm  merged  to 
form  a  single  junction  allowing  for  the  unique  opportunity  to  probe  the  voltage 
inside  of  the  barrier.  It  is  also  strong  evidence  that  the  barrier  in  ion  damage 
junctions  is  much  larger  than  the  slit  size  in  the  implantation  mask. 
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4.3  Series  Arrays  of  Junctions 

The  focus  of  this  work  was  to  fabricate  multi-junction  series  arrays  with 
similar  IcRn  products  for  use  in  ac  voltage  standards.  Arrays  were  fabricated  with 
line  spacing  that  ranged  between  125  nm  -  2  /mi.  The  slits  used  to  define  the 
junctions  ranged  from  30-110  nm  wide.  Arrays  of  4,  10,  20,  50,  and  100  junctions 
were  attempted.  Figure  4.18  shows  a  section  of  a  typical  implantation  mask.  This 
one  in  particular  was  for  a  100  junction  array  with  line  spacing  of  231.7  nm  and 
slit  width  of  50  nm. 

In  an  attempt  to  enhance  uniformity  a  technique  was  developed  to  reduce 
the  effective  width  of  the  bridges  using  ion  damage.  Rectangles  on  the  bridge  edges 
were  written  in  the  PMMA  and  etched  along  with  the  slits  defining  the  junctions, 
so  that  the  edges  would  be  damaged  during  the  implantation.  Figure  4.19  shows 


Figure  4.18:  SEM  picture  of  an  etched  tri-layer  mask  used  for  creation  of  a  100  junction  array.  The  slits 
are  50  nm  wide  with  a  spacing  of  231.7  nm. 
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Figure  4.19:  (a)Top  view  and  47°  tilt  of  a  tri-layer  mask  patterned  so  that  ion  damage  can  reduce  the 
bridges  electrical  width  from  4  to  2  [i m. 


the  top  view  of  the  stencil  used  to  damage  the  junctions  and  bridges  for  a  sample 
of  10  junction  arrays.  Samples  were  made  with  bridges  as  narrow  as  0.5  fj, m. 

Figure  4.20  shows  the  I—V  characteristics  for  ten  junction  arrays  from  3 
different  samples.  Samples  A  and  B  were  fabricated  from  200  nm  YBa2Cu307_,5  on 
LSAT,  implanted  with  2xl013  Ne+/cm2  at  200  keV.  Sample  C  was  fabricated  from 
200  nm  YBa2Cu307_,5  on  AI2O3,  implanted  with  0.5xl013  Ne+/cm2  at  200  keV, 
furthermore  the  bridge  width  was  damaged  down  to  2  n m.  All  of  the  junctions 


Figure  4.20:  I—V  characteristics  for  3  different  ten  junction  arrays  measured  at  78.0  K,  77.3  K  and  84.0 
K  for  A,  B,  and  C  respectively. 


69 


had  similar  slit  sizes  80,  85,  and  90  nm  for  A,  B,  and  C  respectively.  The  inter¬ 
junction  spacing  was  400,  2000,  and  1000  nm  for  A,  B,  and  C  respectively.  The 
I—V  characteristics  are  RSJ-like  and  all  exhibit  flat  giant  Shapiro  steps  at  10  times 
the  voltage  of  a  single  junction.  The  resistance  of  sample  C  is  higher  than  that  of 
A  and  B  as  expected  from  having  the  bridge  reduced  to  2  pm. 

Figure  4.21  shows  the  I  —  V  characteristics  of  sample  C  at  4  different 
temperatures  fit  using  the  RSJ  model  modified  for  thermal  rounding  and  current 
dependant  resistance  as  described  in  the  previous  section.  The  noise  temperature 
calculated  for  the  smallest  Ic  is  89  K  which  is  only  5°  higher  than  the  actual 
measurement  temperature  and  in  good  agreement  with  our  single  junction  results. 

Junction  arrays  were  tested  for  uniformity  using  the  ac  Josephson  effect. 
Appearance  of  flat  giant  Shapiro  steps  is  an  indication  that  the  parameter  spread 
of  the  junctions  in  the  array  is  less  than  10%  [51-53].  Figure  4.22(a)  shows  I-V 
characteristics  for  sample  B,  with  and  without  15  GHz  microwave  radiation.  Giant 
Shapiro  steps  are  visible  at  0.31  mV,  10  times  the  voltage  expected  for  a  single 


Voltage  (pV) 


Figure  4.21:  I—V  characteristics  for  4  temperatures  for  sample  C,  fit  with  a  modified  RSJ  model. 
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Figure  4.22:  Ten  junction  array  I-V  characteristics  measured  at  77.4  K,  with  and  without  15  GHz 
microwave  radiation  (a).  Shapiro  steps  are  visible  at  0.31  mV,  10  times  the  voltage  expected  for  a  single 
junction  at  this  frequency.  Differential  resistance  (b)  measurements  going  to  zero  verify  that  the  steps 
are  vertical,  implying  all  junctions  are  locked  to  the  15  GHz  signal. 


junction  at  this  frequency.  Differential  resistance  was  also  measured  to  verify  that 
the  step  was  the  product  of  all  of  the  junctions.  If  one  junction  has  different 
parameters  it  will  produce  its  step  at  a  different  bias  which  will  have  a  rounding 
effect  on  the  giant  Shapiro  step.  This  is  unacceptable  for  voltage  standards  and 
waveform  generators  which  requires  the  position  of  the  giant  Shapiro  step  to  be 
exact. 

The  microwave  power  dependance  of  the  critical  current  was  also  mea¬ 
sured  for  sample  B  and  is  shown  in  Figure  4.23.  It  shows  a  bessel-like  dependance 
similar  to  single  junctions  however  at  the  higher  order  oscillations  non-uniform 
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Figure  4.23:  Dependance  of  the  critical  current  for  a  ten  junction  array  on  microwave  power. 

critical  currents  of  the  individual  junctions  split  the  peaks. 

Resistance  as  a  function  of  temperature  of  the  junctions  was  measured  to 
determine  the  reduction  in  Tc  from  that  of  the  bulk  material.  Figure  4.24  shows 
measurements  for  Samples  A  and  C.  after  the  bulk  electrodes  go  superconducting 
~  85  K.  The  reduction  in  Tc  is  6.5  K  for  the  higher  implantation  dose  in  comparison 
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Figure  4.24:  Resistance  vs.  temperature  for  samples  A  and  C.  The  electrodes  in  both  films  superconduct 
at  87  K.  The  reduction  in  Tc  is  6.5  K  for  the  higher  dose(sample  A)  and  1  K  for  the  lower  dose  (sample 
C).  The  inset  zooms  in  on  the  transition  for  the  lower  dose  sample. 
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to  1  K  for  the  lower.  The  height  of  the  foot  is  higher  about  twice  as  high  for  sample 
C  which  is  expected  because  the  bridge  width  was  reduced  from  4  to  2  /xm.  The 
transition  of  sample  C  is  smooth  which  indicates  that  Tc  reduction  of  the  material 
in  the  array  is  uniform,  unlike  sample  A  which  shows  transitions  at  both  84.5  and 
80.5  K. 

The  modulation  of  the  critical  current  in  magnetic  field  was  measured  for 
all  three  samples  and  is  presented  in  Figure  4.25  exhibiting  quasi-like  Fraunhofer 
patterns.  Only  half  of  the  pattern  was  measured  for  sample  C  because  the  smaller 
bridge  width  required  more  field  then  the  magnet  could  produce  with  the  power 
supply  used.  The  modulation  period  is  not  constant  most  likely  due  to  unequal 
junction  areas.  In  a  non-uniform  array  the  measured  Ic  is  that  of  the  junction(s) 
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Magnetic  Field  (Gauss) 

Figure  4.25:  Critical  current  modulation  in  magnetic  field  for  samples  A  (a)  78.0  K,  B  (b)  76.0  K,  and 
C  (c)  82.9  K. 
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with  the  smallest  critical  current.  In  the  presence  of  a  field  the  smallest  critical 
current  will  be  determined  by  the  number  of  flux  quanta  penetrating  the  junction’s 
area  as  well  as  it’s  zero  field  value.  The  junction  with  the  largest  area  has  the 
shortest  oscillation  period  as  well  as  the  smallest  critical  current  therefore  the  first 
minima  is  solely  due  to  the  lowest  7c(s)  in  the  array.  At  higher  order  oscillations 
interferences  will  occur  as  the  higher  critical  current  junctions  reach  their  minima. 
The  interferences  are  observed  in  the  data  with  the  exception  of  sample  C  which 
appears  to  be  very  uniform.  Different  levels  of  flux  trapped  in  between  the  various 
barriers  can  also  play  a  big  role  in  determining  the  shape  of  these  patterns.  Table 
4.1  shows  the  average  oscillation  period  of  magnetic  field  oscillation  for  the  three 
arrays  of  Figure  4.25  and  three  single  junctions.  The  data  for  the  single  junctions 
includes  one  fabricated  on  sample  B,  the  junction  from  Figure  4.9  and  a  grain 
boundary  junction  fabricated  from  the  same  material.  These  are  compared  to  the 
three  theories  discussed  previously.  It  is  clear  that  the  ion  damage  junctions  exhibit 
classical  junction  behavior  whereas  the  grain  boundary  junction  closely  follows  the 
planar  junction  models. 


Table  4.1:  Comparison  of  Afi(Gauss) 


Sample 

A B  Actual 

A B  Classical 

A B  Rosenthal 

AjB  Humphries 

Array  A 

17.1 

15.2 

9.5 

4.7 

Array  B 

8.9 

9.7 

3.1 

3.2 

Array  C 

15.0 

15.0 

9.5 

4.7 

Single  B 

9.9 

10.4 

4.2 

3.5 

Single  Fig.  4.9 

14.9 

15.6 

9.5 

4.7 

GB  Junction 

3.0 

8.5 

3.1 

2.8 

Arrays  containing  more  than  ten  junctions  were  not  of  sufficient  unifor¬ 
mity  to  produce  RSJ-like  I  —  V  characteristics  or  flat  giant  Shapiro  steps.  Figure 
4.26  shows  the  I  —  V  characteristics  for  a  50  junction  array  with  200  nm  inter¬ 
junction  spacing.  Large  amplitude  steps  are  visible  at  50  nhf/2e,  however,  the 
differential  resistance  doesn’t  not  reach  zero.  This  is  attributed  to  a  variation  of 
critical  currents  of  the  individual  junctions  within  the  array. 
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Figure  4.26:  (a)  I—V  characteristic  for  a  50  junction  array  measured  at  71.8  K  with  and  without  20  GHz 
microwave  radiation,  (b)  Differential  resistance  measurements  do  not  reach  zero  due  to  non-uniform 
junction  parameters. 


Reproducible  series  arrays  of  Josepshon  junctions  for  ten  or  less  junctions 
exhibit  RSJ-like  I  —  V  characteristics  and  flat  giant  Shapiro  steps  indicating  that 
the  junction  parameter  spreads  are  less  than  10%.  The  critical  current  of  the  arrays 
modulate  in  magentic  field  similar  to  single  junctions  with  a  period  that  can  be 
predicted  with  classical  sandwich-type  junction  theory.  Arrays  of  greater  than  10 
junctions  show  rounded  I—V  characteristics  and  as  well  as  rounded  giant  Shapiro 
steps,  possibly  limited  by  the  EBL  proximity  effect. 


Chapter  5 


MgB2  Josepshson  Junctions 


The  major  focus  of  this  dissertation  was  on  YBa2Cu307_^  junctions,  how¬ 
ever  MgB2  ion  damage  junction  results  by  Kahler  et.  al.  [54]  served  as  an  inspira¬ 
tion  to  extend  this  technology  to  MgB2.  The  initial  results  from  one  of  the  first  few 
samples  was  encouraging  and  exceeded  expectations.  Figure  5.1(a)  shows  I  —  V 
measurements  for  a  single  junction  measured  at  37.2  K,  with  and  without  12  GHz 
microwave  radiation.  The  I  —  V  characteristics  are  RS  J-like  with  an  IcRn  product 
of  75  juV.  Near  Tc  of  the  junction  they  did  not  exhibit  the  subtle  rounding  near  Ic 
as  observed  in  YBa2Cu307--j  junctions.  This  is  an  indication  that  the  rounding  in 
YBa2Cu307_i  junctions  is  indeed  thermal  rounding  and  less  pronounced  at  lower 
temperature. 

Furthermore  at  lower  temperatures  the  MgB2  junctions  did  not  show  a' 
cross-over  to  flux-flow  as  seen  in  the  YBa2Cus07_5  junctions  despite  having  a 
larger  w/Xj  ratio.  For  the  junction  of  Figure  5.1,  w/Xj  ~  8.  This  fits  well  in 
the  framework  of  the  Likharev  SS'S  junction  model.  MgB2  has  a  coherence  length 
that  is  2-3  times  that  of  YBa2Cu307_^,  which  implies  that  L /£  <  3.5  for  a  much 
wider  range  of  temperatures.  Under  microwave  radiation  Shapiro  steps  are  visible 
in  the  I  —  V  characteristic  at  the  expected  voltages:  V  =  nhf  / 2e  and  half  integer 
Shapiro  steps  appear  when  the  junctions  are  measured  at  temperatures  where  Ic  is 
greater  than  0.75  mA.  This  may  be  due  to  magnetic  flux  generated  from  the  bias  or 
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Figure  5.1:  (a)  I  -V  Characteristics  for  a  single  junction  measured  at  37.2  K,  with  and  without  12  GHz 
microwave  radiation,  (b)  Microwave  power  dependance  of  the  critical  current  and  lst-order  Shapiro  step. 

from  a  non-sinusoidal  Josephson  current  relation.  Figure  5.1(b)  shows  the  power 
dependencies  of  the  critical  current  and  1st  order  Shapiro  step.  The  Bessel-like 
dependance  is  a  clear  indication  of  the  AC  Josephson  effect. 

The  critical  current  as  a  function  of  temperature  was  determined  for  a 
single  junction  by  recording  I  —  V  characteristics  at  different  temperatures  (Figure 
5.2).  Like  the  YBa2Cu307_a  junctions  these  junctions  have  an  Ic  cc  (T  —  Tc )3 
near  Tc  and  Ic  oc  (T  —  Tc)2  at  lower  temperatures.  Resistance  as  a  function  of 
temperature  measured  using  a  10  /xA  bias  (Figure  5.2(a))  shows  the  Tc  of  the  bulk 
material  to  be  38.8  K  and  the  Tc  of  the  weak  link  to  be  38.2  K.  This  reduction  in  Tc 
(«  0.6AT)  is  approximately  one  order  of  magnitude  smaller  than  for  YBa2Cu307_5 
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Figure  5.2:  (a)  Single  junction  critical  current(circles)and  resistance  (triangles)  temperature  dependen¬ 
cies.  The  dashed  and  solid  lines  are  fits  with  Ic  oc  (T  -  Tc)2  and  Ic  oc  (T  -  Tc)3  respectively,  (b)  The 
same  data  shown  for  a  temperature  range  near  Tc 

junctions^  6 K)  fabricated  with  identical  conditions.  This  indicates  that  the  Tc 
of  YBa2Cu307_5  is  more  sensitive  to  200  keV  Ne+  ion  damage  than  MgB2  . 

Figure  5.3  shows  the  I  —  V  characteristics  for  a  20  junction  series  array 
measured  at  37.5  K,  with  and  without  12  GHz  microwave  radiation.  A  flat  giant 
Shapiro  step  is  visible  at  20  times  the  value  of  a  single  junction  at  the  same 
frequency.  This  leads  us  to  believe  that  the  spread  in  IcRn  is  small,  however 
rounding  near  Ic  suggests  some  spread  in  Ic.  Differential  resistance  was  measured 
under  the  same  conditions(Figure  5.3(b)).  dV/dl  reaching  zero  confirms  that  the 
step  is  flat  (e.g.  All  junctions  are  locked  to  the  12  GHz  drive  signal).  These  20 
junction  results  are  quite  encouraging  considering  that  this  is  the  largest  array 


Figure  5.3:  (a )  I  —  V  Characteristics  for  a  20  junction  array  measured  at  37.5  K,  with  and  without  12 
GHz  microwave  radiation,  (b)  Differential  resistance  vs.  voltage  measured  under  the  same  radiation. 

size  attempted  and  one  of  the  first  4  samples  fabricated  using  MgB2  .  It  could  be 
a  combination  of  longer  coherence  length  and  more  uniform  material  that  makes 
junctions  from  MgB2  easier  to  fabricate  opposed  to  those  of  YBa2Cu307_<5. 

MgB2  junctions  and  junction  arrays  were  fabricated  that  had  a  similar 
critical  current  dependance  as  those  made  from  YBa2Cu307_5,  suggesting  a  prox¬ 
imity  effect  SN’S  coupling  mechanism.  The  lack  of  observing  a  cross-over  to  flux 
flow  I  —  V  characteristics  in  these  junctions  is  likely  due  to  the  larger  coherence 
length.  The  quality  and  properties  of  MgB2  make  uniform  junctions  and  arrays 
easier  to  fabricate  then  YBa2Cu307_,5. 


Chapter  6 


DC  Washer  SQUID 


One  of  the  problems  with  fabricating  quality  dc  SQUIDS  is  increasing 
the  Resistance  so  that  the  device  can  be  matched  to  other  electronics.  An  attempt 
to  remedy  this  was  to  use  a  junction  array  in  place  of  the  single  junctions  in 
conventional  SQUIDs.  This  was  demonstrating  using  the  mask  shown  in  Figure 
6.1.  10  slits  are  written  in  the  implant  mask  across  each  arm  of  the  SQUID 
resulting  in  2,  10-junction  series  arrays  connected  in  parallel.  Figure  6.2(a)  shows 
the  I -V  characteristic  for  a  NdBa2Cu307-(5  SQUID.  The  curve  is  quite  rounded 
even  near  Tc  ,  which  is  attributed  to  critical  current  non-uniformity.  It  was  later 
found  that  the  material  used  for  this  device  was  of  poor  quality.  Despite  of  the 
non-uniformity  this  device  still  functioned  as  a  SQUID.  A  voltage  modulation  of  10 
fiV  was  observed  at  a  dc  bias  of  115  fiA  (Figure  6.2(b)).  The  period  of  modulation 
is  6.5  mGauss  in  comparison  to  the  value  of  5.5  mGauss  calculated  from  the  sample 
geometry.  The  1  mGauss  discrepancy  is  well  within  the  error  of  the  calibration  of 
the  magnet  used. 

This  demonstrates  a  novel  technique  in  using  arrays  to  increase  the  resis¬ 
tance  for  improved  matching  of  SQUIDs  to  other  components. 
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Figure  6.1:  Mask  used  to  fabricate  DC  washer  SQUIDs.  The  zoomed  region  shows  the  two  arms  of 
the  SQUID  the  lines  crossing  them  represent  the  EBL  defined  slits  in  the  tri-layer  mask  to  define  the 
junctions.  The  square  in  the  hole  of  the  washer  is  an  alignment  mark  which  is  exposed  during  EBL  and 
damaged  with  the  junctions  so  that  it  does  not  superconduct. 


Chapter  7 


Conclusion 


In  conclusion  Josephson  junctions  fabricated  from  electron  beam  litho¬ 
graphy  and  ion  damage  were  fabricated  to  study  their  physics  of  operation  and 
to  improve  uniformity.  The  uniqueness  of  these  junctions  in  comparison  to  other 
types  is  that  they  are  fabricated  in  the  plane  of  a  single  film  with  no  interfaces 
between  different  materials.  This  and  that  the  barrier  is  a  section  of  reduced  Tc 
superconductor  makes  the  physics  of  these  junction  unique.  Future  experiments 
should  explore  the  effect  the  barrier  length  has  on  the  flux-flow  cross  over  as  well 
as  study  the  Fraunhofer  patterns  at  low  temperatures  to  look  for  evidence  of  a 
non-sinusoidal  current-phase  relation. 

Series  arrays  of  10  junctions  for  YBa2Cu307_<5  and  20  for  MgB2  had  uni¬ 
form  junction  parameters  and  exhibited  fiat  Shapiro  steps  which  makes  them  suit¬ 
able  for  waveform  generators,  oscillators  and  voltage  standards.  Further  process 
refinements,  namely  EBL  and  the  incorporation  of  the  junctions  into  a  co-planar 
waveguide  could  substantially  increase  the  number  of  uniform  junctions  that  can 
be  produced. 

Future  work  will  be  to  incorporate  these  junctions  into  a  SQIF  antenna 
device.  This  device  is  less  sensitive  to  parameter  spreads.  Successful  SQIFS  have 
been  fabricated  with  parameter  spreads  as  large  as  30%.  However  SQIFs  require 
junctions  with  large  IcRn  products  of  the  order  of  hundreds  of  microvolts.  This  will 
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be  challenging  but  could  be  possible  by  increasing  the  film  thickness,  decreasing 
the  slit  size  in  the  mask,  increasing  the  damage,  and  optimizing  the  anneal. 

As  the  size  and  price  of  cryo-coolers  gets  smaller  in  the  years  to  come 
there  will  be  a  growing  demand  for  an  HTS  Josephson  technology  from  the  de¬ 
fense,  communications  and  medical  industries.  Although  this  technology  currently 
can  not  provide  large  IcRn  products  at  this  time,  it  has  many  advantages  over 
competing  technologies.  With  a  few  innovations  in  materials  and/or  lithography, 
ion  damage  Josephson  junctions  may  be  able  to  meet  this  demand. 
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